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MHOST VERSION 4.2 USERS MANUAL 
PREFACE 

f x ™** 1 describes the user options available for running the mhost 

eJ f nentS ****9 with beS and shell stxS^uSl 

eouation library ^ Ma ? y °P tlons are available in the constitutive 

ar y> the solution algorithms and the analysis rapabi l i ties . 

An overview of the solution algorithms is presented in Sent ion a * 
Stei 2r C £ 


aiwn^ ly hig ? level of flexibility is achieved in the MHOST code bv 
Section F. 1 explanation of user subroutines is printed in 
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THE MHOST PROGRAM 


A. INTRODUCTION 

.r ia1 The MH 9 S T ^ ogram is a finite element code designed to perform nonlinear 
at M^ S Anai tUrbir ? engint T *** section components. MHOST has been developed 
cL^c^ NAS 3 S 2 ^fiq e 7 arCh Corpo 5 atlon under NASA Lewis Research Cem£ 

o^r. ractor to the united Tectaol °^ 

rvrvr^f innovative mixed finite element procedures implemented in the MHOST 

?^w reSUlt “ s ^9nificantly improved accuracy. Thep^edSr^iS^ 

“? Solu i ion of * s y stQm of nonliiZr 

•S2 S S !S r ^ CtU f eS ' 13356(1 u P° n conc epts Of mixed vaSeS^l^nci?Les 
The standard displacement method is also included as an option. P ?1 

01 «“ “* c-tains three types of inelastic 

(i) siitplified plasticity, 

(if) conventional plasticity, and 

(iii) an advanced viscoplasticity law. 

subr^^r Ur Th^ f n^SaT d rn f terial anisotr °Pi r can be obtained through user 
suorourines . The program contains routines for eigenvalue extraction in 

^aSaS?e a ? modal . anal yses. A single step tij integrator ^tiS 2 
vailable for use in nonlinear transient calculations. Additional 

130 obtaine< ? by adding user subroutines. All the elements 
ailable in the current version are numerically integrated linear or rm»H 
ratic isoparametric elements. All quantities, sSLsS^i st^ S2‘ 

input or internally generated, are represented by nodal points values. ' 

whatT^l2 r to in tw f f e f ° r *** "F* Cale 91x1 its processing is scme- 
v*at smilar to that for a conventional finite element package such as the 

™ -i2r!f al Pur P° se Finite Element code. A few addiScn^SaSStii 22 
required here to control the iterative procedure, as well as the numerical 

Vari ° US VeCt ° rS aSSCCiated «ith the nodal sSessTst^, 

___ if. n thls docume nt, data preparation and interpretation of computed 
results are discussed. A brief note on the underlyiS theo^ and 23222 
mmg is also included. For further theoretical details users should n^for- 

^s^s^a? 601 ® 11031 a " d the Ustln9 of the code ln th ® “HOST 

g^n C S“^eTi. a " alySlS CapabiliUes and i"P“t/output 
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TABLE A. 1 MHOST Analysis Capability 


Element 

Definition 

Options 

Beam 

Plane 

Stress 

Plane 

Strain 

Axi- 

symmetric 

Solid 

Three- 

Dimensional 

Solid 

Three- 

Dimensional 

Shell 

Linear 

Isotropic 

Elasticity 

X 

X 

X 

X 

X 

X 

*1 

Anisotropic 

Elasticity 


X 

X 

X 

X 


Cenposite 

Laminate 






X 

Simplified 

Plasticity 


X 

X 

X 

X 

X 

Elasto- 

Plasticity 


X 

X 

X 

X 

X 

Unified 

Creep-Plasticity 


X 

X 

X 

X 

X 

Stress 

Stiffening 

X 

X 

X 

X 

X 

X 

Centrifugal 

Mass 

X 

X 

X 

X 

X 

X 

*2 

Thermal 

Strain 

X 

X 

X 

X 

X 

X 

Creep 

Strain 

X 

X 

X 

X 

X 

X 

NOTES: 








*1 

Applicable only for linear elasticity. 

*2 

Not separately applicable to the unified creep-plasticity option 
in which these quantities are ccmbined together in the model. 
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Analysis 

Module 

Option 


Quasi- 3 

static 

Analysis 

Buckling 3 
Analysis 

Modal* 3 

Analysis 

Modal* 3 

Superposition 

Transient 

Dynamics 


TABLE A. 1 MHOST Analysis Capability (Cont'd) 


Axi- 

Plane Plane symnetric 
Beam Stress Strain Solid 


Three- Three- 

Dimensional Dimensional 
Solid Shell 


x 


XXX 


X 


X 


X 


XXX 


X 


X 


X 


X x X 


X 


X 


X 


X x X 


X 


X 


X 


XXX 


X 


X 


I^rge Displacement 
Analysis x 


x 


X 


Notes: 


*3 


The effect 
included. 


of stress stiffening and/or centrifugal 


mass nay be 
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TABLE A. 1 MHOST Analysis Capability (Gont'd) 


Iterative 

Solution Quasi-static Transient 

Option Profile Solution Frontal Solution Profile Solution Frontal 

Solution 


Modified*^ 

Newton x x 

Newton 

Raphson 

Update x x 

*5 

Cbnjugate 

Gradient x x 

BFGS 

Update x x 

Secant 

Newton 

update x x 

Line 

Search x x 


x 


x 


X 


X 


X 


X 


Notes: 


*4 

Updating of the tangent array can be stopped at an arbitrary step 
of the Newton-Raphson stiffness update process 

*5 

The line search option is automatically turned on. 
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The MHOST program is written in FORTRAN, with the ANSI-68 standard 
followed as closely as possible. All floating point calculations are 
carried out in double-precision (64 bits/word) on the 32 bits/real-word 
machines (e.g., PRIME, VAX and IBM), or in single-precision (60-64 
bits /word) on the higher precision machines (such as CDC and CRAY) . 

nCREMENEAL-ITERATIVE SOLUTION PROCEDURE 


This section provides a brief discussion of the numerical procedure used 
in the MHOST program. It is not our intention here to discuss the procedure 
in ctetail, but rather to provide a u ser with sufficient background to enable 
him/her to use the MHOST program efficiently for inelastic analyses. A 
cai prehens i ve discussion of the solution procedures used in the MHOST code 
is contained in the MHDST Theoretical ltonual. 

The iterative finite element solution procedure is illustra ted by a flow 
diagram as shown in Figure A.l. Matrices K and F are the conventional 
finite element stiffness matrix and load vector, respectively. Other con- 
ventional finite element matrices appearing in this flow diagram are B, D, 
and Q, respectively the strain-displacement matrix, the stress-strain matrix 
and the diagonalized inner product of shape functions. All quantities 
involved in the computation are evaluated as nodal vectors . Values for 
stresses, strains and displacements are obtained at nodes by the nodal 
strain recovery calculation indicated in Box 3 and the residual load correc- 
tion in Boxes 2, 4 and 5 of the flew diagram. Values at integration points 
are calculated and printed only when requested and these results can be 
compared with those of the conventional finite element method. 

Iterative solutions are utilized even when solving for purely elastic 
response . An iterative solution of elastic response, compared to the stan- 
dard displacement method, produces improved stress results for the same 
accuracy in displacement field. In path-dependent inelastic calculations, 
such mproved response takes on larger importance. 

The user can try various options to conduct an iterative calculation 
procedure. Options include the different criteria to terminate the itera- 
tive calculation as well as numerical quadrature rules to form various 
matrices and vectors. Input data related to the iterative solution are 
documented in Section C.10, for the numerical integration, and D.6 for 
termination of the process. Various additional options have been added to 
the code to improve computational efficiency. It is generally recaimended 
that a user test the options to be able to choose the most efficient com- 
bination for his/her particular problem. 

The iterative solution procedure is embedded in the conventional in- 
cremental finite element control structure as shown in Figure A2. In the 
incremental calculation, the tangent stiffness matrix is formed at each 
increment (Full-Newton) . Within each increment, the iterative algorithm 
performs in the same manner as the Newton- Raphson procedure for solving a 
system of nonlinear equations. The Full-Newton procedure is the default in 
MHOST with the modified-, quasi- and sec ant -Newton iterations as options. A 
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self- adaptive load incrementation option based upon the arc-length method 
can be invoked for computation of the tangent stiffness . When an optional 
method is invoked, the operation in Box 1 of Figure 2 is skipped after the 
specified number of iterations. 

Section E of this document describes the input data associated with the 
incremental procedure used for nonlinear and transient analyses. 

MHOST provides the default values for parameters associated with the 
integration and the iterations, however the user is encouraged to overwrite 
the option defaults to explore the possibilities of the novel computational 
technology implemented into the MHOST code. 
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FIGURE A. 1 Flow Diagram for Iterative Solution Procedure 
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FIGURE A. 2 Nonlinear Finite Element Solution Control Structure 
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B. 

B.l 


DATA INPUT (General) 

Input Structure 

The data input to the program MHDST consists of three blocks . 

1. The PARAMETER DATA block serves to specify and control MHOST 
execution options (e.g., modal analysis, stress stiffening, 
etc.). Also, in this block, upper bounds are specified for 

inode f. p ^ rai ^ ters ( e -9-/ the number of elements, nodes, 
etc. ) to permit the dynamic memory array allocations to be set. 

2 * J he ^ ni J? element model is defined with the MDDEL DATA block. 
In this block of data finite element topology, nodal coor- 
dinates, boundary conditions and loads are specified. 

3. The UO^ffiNTAL DATA block specifies the loading history as well 

t° ^tial and boundary conditions. MHOST has 
no limit on the number of incremental data blocks. 

i n ^ al COndi 1 ti °r S ^ ^ ieriodic loading conditions for the tran- 

?^!Lr? yn ^ niC ca l cu l atl °ns are specified as part of the first 
INCREMENTAL DATA block. F ne rirst 

^^ Q bl > °P k t S t ^ inated with a line containing the keyword END. All 
blocks will contain a number of SECHENTS. Eachsegrrent con- 
sists of a line with a KEYWORD, and if necessary followed bv a 

° £ “if S““- Se9TOnts «“ «te^ in 

once in The Ea ? segment t yP e ma Y appear more than 

h f block. During the second and subsequent ap- 

pearances of a segment in a block, the data entered is processed 
without taking into account previously specified data input, in 
many cases this will mean that the pre^ioSly specifiers is 

SSed earlieT CaSeS ' “ iS alSO P 08 ^ 1 * t0 co * pleme » t data 

ijip £ P f 5Cessor onl Y requires the first four characters of 
wTSS actS 3117 can te < **“ ia ** by its first 

Ton 916 s P ecified dynamic transient carputa- 
m ° n ' t bl ?T k of data segments must be added to the data deck as a 

incremental data block. A special key word to 

li f «r^T-^ l0< ; k 1S u re ^ ired ‘ History dependent loading data 

at ? ach txme increment unless otherwise given by 
means of user subroutines. y * 

The MHOST input data processor reads all of the input data character- 
by-character using FORTRAN A1 format. The input data file is 
processed completely with checks made for possible typing errors and 
data inconsistency, when such a mistake is detected, the program 
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prints comprehensible error messages and terminates execution after 
processing the remaining input data. 

B.2 Line Formats 

In both keywords and numerical data input, items are separated from 
each other by one or more blank spaces. There is no limitation as to 
where data may start or end on a line. The end of a line is also 
considered a separator; hence, every data item must be included in 
its entirety on one line. 

MHOST accepts four types of input lines: 

a . KEYWORD lines 

These lines always start with an asterisk (*), followed by the 
keyword. The star may be preceded by blank, and blanks may also 
be included between the start find the keyword. The keyword may 
be followed by several numbers, separated frcm the keyword and 
each other by blanks. A keyword is recognized by its first four 
characters . 

b. COMMENT lines 

Garment lines may be included by the user to explain various 
parts of his input. MHOST ccmnent lines start with a C, 
optionally preceded and/or followed by blanks. The text on a 
ccmnent line is reflected literally in the output file. 

Otherwise ccnment lines have no influence. Garpletely blank 
lines are treated as carment lines. 

c. DATA lines 

These lines contain only numbers. DATA lines cannot start with 
a * or a C. Numbers are separated from each other by blanks. 

d. The TITLE line 

The first line of the input file is considered to represent the 
title of the problem. This line is printed at the top of every 
page. 

B.3 Numerical Data 

Numerical input data consists of integer and real data. 

An INTEGER consists of a number of digits, optionally preceded by a 
plus ( + ) or a minus (-) sign. Blanks are not allowed within the 
number. 
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REAL can be specified in various ways. Hie simplest way is to use 
the same format as used for integer specification. Alternatively, a 
REAL can be specified with a decimal point. The decimal point can be 
positioned anywhere within the range of digits, or before or after 
the complete series of digits, but never before the plus (+) or minus 
(-) sign, if present. v ' 

If the REAL consists solely of a decimal point (optionally preceded 
by a sign), the REAL has the value zero. This means that ' . ' »+» 

' ~ ' • all will be read as zero . ' 

Finally, a REAL can be specified with an exponent. The exponent 
consists of the letter E followed by an integer, and is attached at 
the end of a real as specified above. If only the exponent is 
specified, it is assumed that it is preceded by a real of value 1 
(i.e. E - 10 = 1.0E-10). If a line contains more numeric entries 
• reqU ^ red ' ^ rernainin g numeric entries will not be read in. In 
addition, it is possible to exclude the E in the exponent if the 

is preceded with either a + or - (i.e., 2.3 +6 
^.JEb). If a line contains fewer numeric entries than required 
the remaining numeric entries will be set equal to zero. 

B.4 Continuation T.in^g 

° f 1 i ntegerS ^ reals on a line is tied to certain maxima: 

t-p that 16 ■ Lnte 3 ers or 8 rea l s can be entered on a single line. 

i + .Z? th integers 91x1 rBals appear on a line, the limit can be forrnu- 
latea as z 


(NO. 


of integers ) *5 + (No. of reals)*10 < 80. 


If T JJ F Ut u r ? COrd exceeds the above limit, continuation lines are 
used, which have the same specification as regular numerical data 
lnies. Continuation lines are only used following numerical data 
lines in which no space for further data is provided. 


B.5 Files 


The MHOST program uses six files for input and output. The files and 
their contents are listed below: 


Unit 


Contents 


1 

5 

6 

8 

11 

19 


Log file 
Input data 

Input Echo, Error messages, Iteration 
data, Processing Information and Results. 
Restart data, (binary form) 

Scratch file. 

Post-processing file 
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The names of the above files are system dependent. 


B.6 Progr aiming Considerations 

The parameter data lines in the input file are read by SUBROUTINE 
DATIN1. SUBROUTINE DATIN1 sets up flags and parameters used inter- 
nally to control execution. 

MHOST utilizes dynamic core allocation for memory management. Core 
allocation is carried out in SUBROUTINE INITI1 which is entered frcm 
the model data reader SUBROUTINE BULKIN. SUBROUTINE INITI1 gives 
values to pointers for the mesh data storage. Subsequently, 
SUBROUTINE DATIN2 is entered to read the model data. Individual 
data input modules are executed with respect to the key wards. 

Most trivial input data errors are detected by the code itself. 
MHOST also performs a cross-reference check of the parameter and 
model data for consistency. 

It is advisable for a user to go through the input data whenever 
something goes wrong. For the post-mortem, the MHOST Systems Manual 
should be consulted. 

Once the input of the model data has been completed, one of the 
finite element driver routines is entered. These routines are 
listed in Figure B.l. The coding strategy is that a new analysis 
capability is added to the system by introducing a new driver 
routine. Each driver routine uses common finite element and con- 
stitutive equation libraries. 

The incremental data input is processed by SUBROUTINE INCRIN which 
calls SUBROUTINES DATIN3 and DAT0H1 . In SUBROUTINE DATIN3 , the 
incremental data lines are read and added to the existing in-core 
database. The automatic and adaptive load incrementations are 
processed in this subroutine. SUBROUTINE DAT0H1 generates the 
report on the incremental update of loading data. 

The finite element solution is printed by SUBROUTINE PRINOU which is 
used by are driver routines. 

When requested, the post-file generator SUBROUTINE POSTOU is entered 
and a formatted file is written for graphics post-processing. The 
format is consistent with the MARC general purpose finite element 

* 

package and can be processed by MENTAT with a MHOST enhancement. 


*MENTAT is a proprietary graphics/finite element modeling software 
package available frcm MARC Analysis Research Corporation. 
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The finite element driver subroutines execute the finite elsnent 
analysis modules . The execution is terminated by a driver routine 
only when the iterative process fails to converge within the 
specified number of iterations. The execution time of each opera- 
tion is reported in the line printer file (FORTRAN unit number 21) . 


Major analysis modules are SUBROUTINE ASSEM1 which constructs the 
element stiffness matrices, and assembles the global stiffness 
matrix and load vector. SUBROUTINE ASSEM2 performs frontal (wave- 
front) solution of equations. SUBROUTCCN ASSEM3 is used to generate 
the coefficient matrices for transient time integration. Nodal 
strain and stress arrays are computed in SUBROUTINE ASSEM4 No 
diagnostics are generated in these routines. 


The above finite element analysis modules control the subroutines 
which contain both element and constitutive equation libraries 
SUBROUTINE DERIV controls the element library aTthe elemiit 

f BR0OTIN ^^2P EIj 1S "s* 1 111 ^ element lo °P to construct element 
arrays . SUBROUTINE NDDSTR is called by SUBROUTINE STRESS and super- 
vises the execution of element library routines. SUBROUTINE NODSTR 
controls the loop over the nodal points and, for shell elements, the 
integration layers in conjunction with the pre-integration. 

The subelement iteration loop is entered from the residual force 

^™ T l0 ° P in SUBR0UTINE ASSEM4. The main control routine is 
SUBROUTINE SUBFEM, which calls the subroutines needed for a caiplete 
mixed iterative solution in a subelement region. Execution timing 
is reported in the subelement solution process . For the subelement 
solution, the element library subroutines are constructed slightly 
different from the global solution. The constitutive equation 

^ ^ S u a ? ces ? ed via SUBROUTINE STRESS in exactly the same manner 
as the global solution. 
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MAIN Program - BLOCK DATA SUBPROGRAM 

SUBROUTINE HOST 

SUBROUTINE DATIN1 

Reads parameter data 

SUBROUTINE BULKIN 

Allocates core, reads model data and echoes the input data 

- SUBROUTINE INCRIN 

Reads the incremental data and echoes the input data. This 
routine is called from the analysis driver routines, inside of 
the incremental solution loop. 

SUBROUTINE STATIC 

Drives the incremental iterative solution for quasi-static 
analysis (with the band matrix option) 

SUBROUTINE DYNAMT 

Drives the direct time integration for transient analysis 
SUBROUTINE MDOAL 

Drives the eigenvalue extraction subsystem for modal analysis of 
linear elastic structures 

- SUBROUTINE SUPER 

Drives modal superposition for linear elastic structures. Mode 
shapes and frequencies are obtained from SUBROUTINE M3DAL. 

SUBROUTINE FREDQM 

Drives frequency domain analysis for linear elastic structures. 
SUBROUTINE FRONTS 

Drives the incranental iterative solution for the quasi-static 
analysis (with the frontal solution option) 

SUBROUTINE BUCKLE 

Drives the linear buckling analysis. 

SUBROUTINE PRINCXJ 

- SUBROUTINE PQSTOU 

Generate reports on the line printer file and the formatted post- 
processing file. These routines are called from the analysis driver 
routines. 


FIGURE B.l MHOST Finite Element Driver Routines 
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C. PARAMETER DATA 


sr5 ld ~ 

where the keyword is used as title for the segment. ^ segment. 

The following 32 keywords can be used in Version 4.2: 


SECTION C INDEX 


KEYWORD 

♦ANISOTROPY 

♦BANDMATRIX 

♦BASEEXC ITATIGN 

♦BEAMSECTION 

*BPGS 

♦BOUNDARY 

♦BUCKLE 

♦CENIRIFUGALMASS 

♦COMPOSITE 

♦CONJUGATEGRADIENT 

♦CONSTITUTIVE 

♦CREEP 

♦DASHPOTS 

♦DEFORMATION 

♦displacemenimeihod 

♦DISTRIBUTEDLQAD 

♦DUPLICATENODE 

♦DYNAMIC 

♦ELEMENTS 

♦EMBED 

*END 

*FINITE_STRAIN 

*PQLLCWER_FORCE 

♦FORCES 

♦FRONTALSOLUTION 

*HARDENING_SIOPE 

♦HARMDNICLQAD 

♦LARGEJDISPLACEMENT 

♦LINESEARCH • 

♦LOUBIGNAC 

♦MASSES 

♦MODAL 

♦NODES 

♦NOECHO 

♦OPTIMIZE 

♦PERIODICLQADING 


SECTION NO. 

PACE 

C.15 

22 

C.35 

26 

C.50 

30 

C.40 

27 

C.43 

28 

C.3 

19 

C.17 

22 

C.53 

30-31 

C.51 

30 

C.54 

31 

C. 19 

23 

C. 14 

22 

C.45 

28-29 

C.37 

26-27 

C.41 

27-28 

C.20 

23 

C.32 

25 

C. 12 

21 

C.l 

17 

C.38 

27 

C.58 

33 

C.61 

33 

C.62 

33 

C.6 

19 

C.36 

26 

C.59 

32-33 

C.49 

29-30 

C.60 

33 

C.47 

29 

C.10 

20 

C.46 

29 

C.16 

22 

C.2 

18 

C.55 

31 

C.13 

21 

C.31 

25 
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♦POST 

C.33 

25-26 

♦POWERSPECTRUM 

C.56 

31-32 

♦PRESSURE 

C.58 

32 

♦PRINTSETS 

C.7 

20 

♦PULSELQAD 

C.57 

32 

♦REPORT 

C.21 

23 

♦RESTART 

C.8 

20 

♦SCHEME 

C.30 

25 

♦SECANINEWION 

C.48 

29 

♦SHIFT 

C.42 

28 

♦SPRINGS 

C.44 

28 

♦STIFFENING 

C.52 

30 

♦STRESS 

C.ll 

21 

♦TANGENT 

C.22 

24 

♦TEMPERATURE 

C.9 

20 

♦THERMAL 

C. 18 

23 

♦TRANSFORMATIONS 

C.5 

19 

♦TYING 

C.4 

19 

♦UBOUN 

C.34 

26 

♦UCOEF 

C.26 

24 

♦UDERIV 

C.29 

25 

♦UFORCE 

C.24 

24 

♦UHOOK 

C.28 

25 

♦UPRESS 

C. 27 

24 

♦UIEMP 

C.25 

24 

♦UIHERM 

C.23 

24 

♦WKSLP 

C.39 

27 
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C.l * ELEMENTS 


This parameter must be present to allocate memory for element con- 
nectivity. At least one element type data line is required. 

WARNING: The execution will be terminated if this data line is not 

present. 

Parameter 1 : Upper bound to the number of elements . 

The user should preferably specify the exact number of elements in 
the mesh as the parameter value. The core allocation and print line 
control are based on this parameter. 

After the keyword line, a continuation line is used to define the 
element type to be used in the analysis. In the current version of 
MHOST finite element models are limited to a single element. One 
integer is to be specified on the following line, designating the 
element type. 

The following element types are currently available: 

3: 4-node plane stress element. Node and integration points are 

numbered counter-clockwise. Analogous to MARC element type 3. 
Five element properties: thickness. Young's modulus, Poisson's 
ratio, thermal expansion and density. 

7: 8-node solid element. Node and integration points are numbered 

counter-clockwise. Analogous to MARC element type 7. Five 
element properties: dimity. Young's modulus, Poisson's ratio, 

thermal expansion and density. 

10: 4-node axisyrmetric solid element. Node and integration points 

are numbered counter-clockwise. First coordinate axis is 
symmetry axis. Analogous to MARC element type 10. Five ele- 
ment properties: dummy. Young's modulus, Poisson's ratio, 

thermal expansion and density. 

11: 4-node plane strain element. Node and integration points are 

numbered counter-clockwise. Analogous to MARC element type 11. 
Five element properties: thickness. Young's modulus, Poisson's 

ratio, thermal expansion and density. 

75: 4-node thick/ thin shell element. Nodes and integration points 

are numbered counter-clockwise. Analogous to MARC element type 
75. Five element properties: thickness. Young's modulus, 

Poisson's ratio, thermal expansion and density. 

98: 2— node Timoshenko beam element. Analogous to MARC element type 
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101: 9-node quadratic Lagrangian plane stress element. Node and 
integration points are numbered counterclockwise with the 
center node defined last. Other definitions are identical to 
those in element type 3. 

102: 9-node quadratic Lagrangian plane strain element. Node and 
integration points are numbered counter-clockwise with the 
center node defined last. Other definitions are identical to 
those in element type 11. 

103: 9-node quadratic Lagrangian axisymmetric element. Node and 
integration points are numbered counter-clockwise with the 
center node defined last. Other definitions are identical to 
those in element type 10. 

151: 4-node plane stress element based on an assumed strain formula- 
tion. Node and integration points are numbered 
counterclockwise. Six entries for element properties include 
the thickness, Young's modulus, Poisson's ratio, thermal expan- 
sion coefficient, mass density and (optional) shear modulus, in 
that order. 

152: 4-node plane strain element based on an assumed strain formula- 
tion. Node and integration points are numbered 

counterclockwise. Six entries for element properties include a 
dummy parameter. Young's modulus, Poisson's ratio, thermal 
expansion coefficient, mass density and (optional) shear 
modulus, in that order. 

153: 4-node axisymmetric element based on an assumed strain formula- 
tion. Node and integration points are numbered 

counterclockwise . Five entries for element properties include 
a dummy parameter, Young's modulus, Poisson's ratio, thermal 
expansion coefficient and mass density, in that order. The 
shear modulus should not be specified separately, since this 
would make little sense for an axisymmetric problem! 

154: 8-node solid element based on an assumed strain formulation. 
Node and integration points are numbered in accordance with the 
convention shown in Figure G.10. Six entries for element 
properties include a dummy parameter. Young's modulus, 
Poisson's ratio, thermal expansion coefficient, mass density 
and (optional) shear modulus, in that order. 

Detailed element descriptions are provided in Section G. 

C.2 *N0DES 

This keyword must be present to allocate memory for the nodal point 

coordinate values. The number of entries at each nodal point is 

calculated in the code. 
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Mints 1 *?? 1 ??? 611 t ? at user specif y the exact number of nodal 
& ln tbL he n^: C ° re all ° catio " and Prt-t l-e control are 

Parameter Is Upper bound to the number of nodes. 

C.3 * BOUNDARY 

Stf *“e»SS! ?nr?? es men ! OI i Y £or n0dal dis Placement constraint 
data. Except fcpr a few special cases, at least the same number of 

boundary conditions as the rigid body modes must be specified. 

Insufficient specification of the boundary conditions can result in 
a singular stiffness matrix. result m 

Parameter Is Upper bound to the number of boundary conditions. 

C.4 *TYING 

r^l d dat a ^ aC ™ iad * 

™ C °X!i train u is ^P 03 ^ if the corresponding data segment is 
provided m the model data section. No diagnostics are generated. 

Parameter 1: Upper bound to the number of tying equations 

Parameter 2: Upper bound to the number of degrees-of-freedan 
involved in a tying equation. 

C.5 *TRANSPORMATICNS 

f ° r coordinate transformations of 
nodal degrees-of-freedan. The angles of rotation are specified bv * 
data segment accompanied by the same keyword in the modefdata. ** 

^^s formations take place only when the data is supplied in the 
model data section. No diagnostics are generated. 

Parameter 1: Upper bound to the number of transformations. 

C.6 * FORCES 

This keyword allocates memory for concentrated nodal forces The 

i!??? to d6£ine the data segment?^ iSETfoSS 

uiput in the model data and incremental data sections. 

*** applied wnless the data segments appear in the 
model data and incremental data sections. 
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Parameter 1: Upper bound to the number of nodal degrees-of-freedan 

having non-zero concentrated forces. 

C.7 *PRINrSETS 

This optional keyword does not need to be specified unless more than 
10 print options are to be chosen, or if the iterative updates for 
incremental displacement and reaction forces are to be printed. 

Parameter 1: Maximum number of sets to be used in print 

specification. Default is 10 sets. 

A negative parameter value invokes the printing of iterative dis- 
placement updates. 

C.8 * RESTART 

If this line is included, a previous analysis is to be continued. No 
other parameter or model input will be used or is necessary in a 
restart input file. Only the *END line must be included. The 
binary restart file must be available on a file assigned to unit 8. 

C.9 ^TEMPERATURE 

No parameters. Indicates that nodal temperatures can be stored and 
used in the calculations . 

C.10 *LOUBIGNAC 

This keyword invokes the mixed finite element solution algorithms 
and selects the integration schemes for stiffness, strain recovery 
and residual terms. With these opitions the type of mixed interpola- 
tion procedure may be controlled. 

The following parameters are effective to the linear two and three- 
dimensional elements in MHOST Version 4.2. 

For reduced integration, strain sampling and numerical quadrature 
are carried out at the centroid of the linear elements. For full 
integration, strain sampling and quadrature are carried out at two 
points in each element coordinate direction. 

For the trapezoidal, rule strain sampling and numerical quadrature 
are carried out at each comer of the quadrilateral element. 

For selective integration the shear strain stiffness component is 
evaluated at the element centroid, with other strain components 
evaluated either by the full or trapezoidal integration rules. 

The trapezoidal rule automatically implies use of selective integra- 
tion in the present version of MHOST. 
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Parameter 1: 1. Reduced integration for strain recovery 

2 • Full integration for strain recovery 

3. (Default) Trapezoidal rule for strain recovery 

Parameter 2: 1. (Default) Full integration for residual force 

calculation 

2. Reduced integration for residual force calculation 
Parameter 3: 1. Full integration for stiffness natrix assembly. 

2. Selective integration for stiffness matrix 
assembly. 

3. (Default) Selective integration with element 
Cartesian coordinate transformation for stiffness 
matrix assembly. 

*IOUBIQ'IAC 3 1 3 is the system default and recommended for most of 
the applications. 

C.ll 'STRESS 

This option is used to allocate the memory for storing the nodal 
constraint for the stress components. It is not reccrtmended to use 
this option in practical computations . When this option is invoked 
erroneous results or* divergence iray occur* 

Parameter 1: Upper bound to the number of nodal stress components 

that are to be prescribed. 

C.12 'DYNAMIC 

This option invokes the direct time integration procedure. 

C.13 'OPTIMIZE 

This options is used to flag the bandwidth optimization procedure. 

Parameter 1: Number of iterations in the Cuthill-MoRee band width 

optimization algorithm. The default value is 10. 

When this option is invoked, the bandwidth is reported on both line 
printer and log files. Note that the bandwidth is calculated based 
on the node numbers and assuming a single degree-of -freedom per 
node. The subsequent analysis is performed using the first smallest 
nodal numbering achieved during the iteration. The actual bandwidth 
reported by the solver printed on the line printer file is the 
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actual number of non-zero entries in the banded/profile equations. 
The relation is: 

(Actual Hal f -bandwidth ) = 

(half-bandwidth reported by the optimizer) 
x (number of degree-of-freedcm per rode) - 1 


C.14 *CFEEP 

This option is invoked if the conventional creep model is being 
used. 

C.15 * ANISOTROPY 

This option is used to flag anisotropic material behavior. 
Anisotropy may be included for elasticity, yield surface, and ther- 
mal strains. Material axes orientation must be defined in the model 
data. 

C.16 *MDDAL 

This option is used to flag the eigenvalue extraction option. 
Currently the subspace iteration method is utilized. 

Parameter 1: Number of modes to be extracted, defaults to one if 

not given. 

Parameter 2: Number of subspace dimensions for the eigenvalue 

extraction. 

(A sufficient dimension is reserved as a default) 

Parameter 3: Increment number at which eigenvalue extraction is to 

be performed. (Default = 0). 

If the results from the modal extraction are to be used for a linear 
dynamics analysis using modal superposition, a continuation line 
must follow, defining 

Parameter 1: Number of dynamic modes to be used for mode 

superpos ition . 

C.17 *BUCKLE 

This option is used to flag the linear elastic buckling analysis . 

Parameter 1: Number of buckling modes to be extracted. 

Parameter 2: Number of subspace dimensions for the eigenvalue 

extraction 

(A sufficient value is assigned as a default) 
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C.18 


C.19 


C.20 


C.21 


Parameter 3: Increment number at which eigenvalue extraction is to 

be performed. (Default = 0). 

♦THERMAL 

This option is used to flag temperature dependent material 
properties. These properties are given via user subroutine UTEMP. 

♦CONSTITUTIVE 


This option selects the constitutive relation to be used. Note that 
the default is von Mises plasticity as indicated below, not linear 
elasticity. For linear elastic analysis, the user needs to flag the 
correct option here so as to avoid unnecessary ccmputations . 


When the *COMPOSITE option is used for shell elements, 
request for the constitutive equations are ignored. 


all the 


Parameter 1: Option Type 


0 . Li ne ar elasticity. 


1. Simplified plasticity (secant elasticity) model. 

2. Conventional plasticity (default). 


3. Walker's creep plasticity model 

When either the secant plasticity or the conventional plasticity 

J*® limit of elastic response must be specified by 
the WORKHARD data segment (see Section D.18). If such data is 
missing, the secant elasticity model will result in a singular stiff- 
ness matrix and the conventional plasticity model will assume an 

astronomical number (10^) as the yield stress. 


When the Walker's creep-plasticity model is flagged, the user must 
provide parameter values through user subroutine WALCCN. An example 
of this subroutine is included in this document. Also the MHDST 

program is shipped with this exanple. See Section F.9 for additional 
information . 


♦DISTRIBUTEDLQAD 

The distributed load option (including gravity, centrifugal and 
traction loading types) are flagged. No parameter option is needed. 

♦REPORT 


This option controls the output of nodal solutions. 
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C.22 


C.23 


C.24 


C.25 


C.26 


C.27 


Parameter Is Print interval 

Default is nodal solutions specified by the ‘PRINTOPTION ( Section 
D 141 to be printed at the erri of every increment. When thus °E*l°n 
?s fla^dltS solution is printed only at the intervals specrfied 
by the positive integer parameter. 

★tangent 

This option invokes the irodified Hewton-Raphson iteration procedure. 

Parameter 1: Iteration cycle at which updating of the tangent 

stiffness matrix is halted in each increment. 

When this option is flagged, the tang™* is 1 

sDecified iteration count is met, after which no further reassemoiy 
of^he stiffness matrix occurs. If this option is not flagged, the 
tangent stiffness matrix is updated at every iteration cycle except 
for purely linear elastic problems. 

★UTHERM 

This option invokes the call to the user subroutine UTHERM. No 
parameter. When this option is flagged, the user subroutine must be 
upcLted and compiled prior to theexecution.This °Pti° n prevents 
the user subroutine from interfering with other problems. 

♦UPORCE 

This option invokes the call to the user subroutine UTORCE. . » 

parameter. The user subroutine must be in the load module whe 
MHOST is executed with this option. 

★UTEMP 

This option invokes the call to the user subroutine UTEMP. No 
parameter. The user subroutine must be in the load module when 
MHOST is executed with this option. 

★UCOEF 

This option triggers the call to the user subroutine UOOEF. No 
parameter. The user subroutine must be in the load module w 
MHOST is executed with this option. 

★UPRESS 

This option flags the user subroutine call to UPRESS. » P^f^l 
The user subroutine must be in the load module when MHOST 
ecuted with this option. 
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C.28 *UHOQK 

This option invokes the call to user subroutine UHOOK. No 
parameter. The user subroutine must be in the load module when 
MHOST is executed with this option. 

C.29 *UDERTV 

This option flags the call to user subroutine UDERIV. The user 
subroutine must be in the load module when MHOST is executed with 
this option. 

C.30 *SCHEME 

The time integration scheme is altered by invoking this option. 
After the keyword, a free format numeric data line is used to 
specify integration parameters. The second data line consists of 
three reals. Those ares 

Real 1 Reserved for future extension (set 0.0) 

Real 2 Parameter b in Newmark family of single step sc hares 

Real 3 Parameter g in Newmark family of single step schemes 

Note that the default value for b and g are 0.25 and 0.5, respec- 
tively, resulting in the average acceleration scheme. For possible 
combination of these parameters, see Theoretical Ma nual 

C . 3 1 *PERIODICLOADING 

This option is used to invoke periodic nodal force loadings and/or 
displacement boundary conditions. TVjo integer parameters. 

Parameter 1: Maximum number of displacement conditions. 

Parameter 2 ; Maximum number of nodal force conditions . 

The loadings and/or boundary conditions must be defined in the 
incremental data. 

C.32 *DUPLICATENODE 

The duplicated node option is invoked by entering this parameter 
At duplicated nodes the displacement degrees— of— freedom are 
degenerated to the master node. Strains and stresses are still 
evaluated separately at each node. 

Parameter 1: Maximum number of the duplicated nodal points. 

C.33 *PQST 
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This option invokes the generation of post-processing file 
(formatted on FORTRAN unit number 19) for plotting. The contents 
and format of the post-processing file are discussed in Section J. 

Parameter 1: Post file generation interval. 

C.34 *UBOUN 

This parameter is used to flag the call to user subroutine UBOUN. 
The user subroutine must be in the load module when MHOST is ex- 
ecuted with this option. 

C.35 *BANDMATRIX 

The symmetric profile matrix solver is invoked by using this option. 
Note that this option is the default. 

No parameter. 

The bandwidth optimization option can be used. No out-of-core 
solution capability is supported in the current version of MH OS T. 

C . 36 *FRQOTALSOLOTIGN 

The optional nonsynmetric fronted, solution subsystem is entered by 
flagging this parameter. This is an out-of-core solution and util- 
izes a direct access device connected to FORTRAN unit number 14. 

No parameter. 

This option cannot be used with the RANDMATRIX option. No wavefront 
optimization option has yet implemented in MHOST. 

This option cannot be invoked when the dynamic analysis or the 
eigenvalue extraction option is flagged. For certain cases when the 
parameter specification is contradictory, the program generates 
appropriate error messages. Otherwise, this option will either be 
ignored or results in execution error s . 

This option is not operational on IBM systems. 

C.37 *DEFORMATIGN 

This option is used to invoke eigenvalue extraction for the dis- 
placement stiffness matrix. Two integer parameters can be 
specified. 

Parameter 1: Load step at which the deformation modes are 

extracted. 

Parameter 2: Number of eigenvalues and eigenvectors to be 
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C.38 


C.39 


C.40 


C.41 


extracted. 


par f Teter is i^ot given, the program automatically extracts 
the longest eigenvalue at the undeformed initial state. 

This option is extremely useful to examine the characteristics and 
quality of elements and meshes. See 0. C. Zienkiewicz and q 
to thTV 198 ^' ”'^ ie Penalty Function Method and its Applications 

Voi ^q 1 )^ 011 ? f Boundar y Value Problems", ASME AMD. 

vol. 51, pp 157-179. The information obtained from t his calculation 
is equivalent to the results of the patch test. calculation 

♦EMBED 


This keyvord flags the execution of the subelement iteration option. 

^ specified to choose the integration 
algorithm for the subelement stress. 

Parameter 1: 0 For area average trapezoidal integration. 

1 For nodal trapezoidal integration. 

♦WKSLP 


This keyword targets the user subroutine WKSLP to supplv stress- 

PTTJ&JZ £S£: plastic *«• 

♦BEAMSECTION 

! US t d t Z SpeC ^ fy 0ne or 100110 sets of cross-sectional 
properties for elastic beam elements * 

Parameter Is Upper bound to the number of beam element cross- 
sectional property sets. Defaults to 1. 

♦DISPLACEMENIMETHOD 

This option disables the mixed iterative solution algorithm. The 
DrabW 100 ! 3 ! 1 dis P^ acement method is invoked for linear elastic 
S^SSed PrOC Ure for ^^tic computations is not 

No parameter. 

JJflJ ****1 K Pt l? n is flagged , the stress at a quadrature point is 
calculated basal on the displacement method and isused to generate 

f S rCG VeCt ° r 111 equilibrium iteration. The total 
d0e l - n0 ^ ^one precisely zero due to round off. 
Therefore, it is advisable to specify a fictitious relative residual 
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tolerance to avoid the ERROR logo print out and discontinuation of 
the execution. 

The strain and stress are recovered at nodes even when this option 
is used for the solution. The MHOST code does not support any 
option to directly recover strains and stresses at quadrature 
points. 

C.42 * SHIFT 

This option is used to invoke the power shift capability, allowing 
the specification of one or more frequency ranges for eigenvalue 
search in linear dynamic analyses. 

Parameter 1: Upper bound to the number of power shifts (frequency 

ranges) specified for eigenvalue search. 

C.43 *EFGS 

The quasi-Newton, inverse BFGS update process is invoked by this 
parameter. One integer parameter is required: 

Parameter Is Number of trial vectors to be stored in-core for the 
BFGS update. The default value is 10. 

The default algorithm for the iterative solution of the mixed system 
of equations is the Newton-Raphson tangent stiffness assembly with 
factorization at each inelastic iteration process. BFGS avoids 
these operations and directly updates the inverse of the tangent 
stiffness. This requires only a few vector multiplications along 
with back substitution. 

This option is potentially far more economical than the default 
procedure. It is advisable to flag this option when the standard 
procedure is unacceptably slew or fails to converge. 

The initial tangent modulus is used in this solution algorithm 
(often referred to as the KT-1 method) . The *TANGENT option is 
ignored when the BFGS parameter is specified. 

This option cannot be used with the secant Newton option. BFGS is 
applicable for quasi-static and direct integration dynamic problems. 

C.44 * SPRINGS 

This option allows additional stiffness terms at a node to be added 
directly into the stiffness matrix. 

Parameter 1: Upper bound to the number of added stiffness terms. 

C.45 *DASHP0TS 
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Used to specify additional viscous damping terms at a node, to be 
added directly into the damping matrix. 

Parameter Is Upper bound to the number of added viscous damping 
terms. 


C.46 *MASSES 

Used to specify additional mass terms at a node, to be a d ded 
directly into the mass matrix. 

Parameter Is Upper bound to the number of added mass terms. 

C.47 *LINESEARCH 

This parameter invokes the line search operation in the iterative 
process . 

No parameter. 

The maximum number of search distance iterations is internally set 
to five. When MHOST fails to detect the zero point, an appropriate 
message is printed and execution continues. 

The current version of the MHOST code does not support this option 
when used in conjunction with the arc-length method. Any other 
solution option can be accelerated by invoking this option. 

The operation involves multiple evaluation of the residual vectors 
and can be potentially time-consuming for some problems. 

C.48 *SECANTOEWTON 

When this option is flagged, secant Newton iteration is invoked. 

No parameter. 

The structure is almost identical to the BPGS update. Hie same 
potential economy is expected by utilizing this procedure. 

This option cannot be used with the BPGS option. SECANINEWrCN is 
available for quasi-static and direct integration dynamic problems . 

C.49 *HARM3NICLQAD 

This option is used to invoke harmonic nodal force loading for 
dynamic analysis. One integer parameter. 

Parameter Is Upper bound to the number of degrees-of-freedcm at 
which harmonic loads are specified. 
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In addition, a continuation line must be provided with a single real 
parameter, specifying: 

Parameter 1 (REAL): Frequency of excitation, in cycles/time. 

C.50 *BASEEXCITATION 

This option is used to invoke harmonic base excitation at one or 
more base locations. One integer parameter. 

Parameter 1: Upper bound to the number of degrees-of-freedcm at 

which harmonic base excitations are specified. 

In addition, a continuation line must be provided with a single real 
parameter, specifying: 

Parameter 1 (REAL) : Frequency of excitation, in cycles /time. 

C.51 *CCMPOSITE 

The composite laminate model is used for defining the characteris- 
tics of shell element type 75, when this parameter data is 
specified. This option overwrites the element parameter data (See 
Table ELIB.l) and, thus, must appear after the element type defini- 
tion (*ELEMENIS) parameter line. No parameter. 

Note that this option requires the nodal material data definition in 
the model data section using the LAMINATE option. 

This option was developed from the formulation of M. D. Minich and 
C. C. Chamis (1975), Analytical Displacements and Vibrations of 
Cantilevered Unsymmetric Fiber Composite Laminates, NASA Technical 
Memorandum NASA IMX-71699. The underlying assumption is that the 
material exhibits linear elastic response. Any requests for inelas- 
tic constitutive laws are ignored in the code when this option is 
invoked. Creep and thermal strain options are available in conjunc- 
tion with this option. 

C.52 * STIFFENING 

When this option is requested, stress stiffening due to initial 
stresses is invoked. 

Parameter: Increment number at which stiffening is to be included. 

All incrsnents following this increment will also have 
stiffening. (Defaults to increment number one). 

C.53 *CENTRIFUGAIMASS 
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When this option is requested centrifugal mass stiffness effects due 

to angular velocity are included. (Default to increment number 
one) . 

Parameter Is Increment number at which centrifugal mass stiffness 
effects are to be included. All increments following 
this increment will also have centrifugal irass effects 
included. 

The angular velocity is defined with the DISTRIBOTEDLOAD model data 
input. Hie axis about which the angular velocity occurs is defined 
with the BODYPORCE model data input. 

C.54 *CCNJUGATEGRADIENT 

This parameter invokes the conjugate gradient method for the mixed 
finite element equations. The stiffness equations for the displace- 
ment method are used as the preconditioner. The line search is 
automatically turned on with this option. Not usable with thequasi- 
Newton and secant-Newton options. The algorithmic detail is found 
in the paper by Nakazawa, Nagtegaal and Zienkiewicz ( 1985), 
"Iterative Methods for Mixed Finite Element Equations ", ASME AMD 


C.55 * NOECHO 


This option suppresses the direct echo print of the bulk data input. 
This option reduces the amount of line printer output for large 
finite element models. 

C.56 *PCWERSPECTRUM 

Used to invoke the driver for frequency-domain analysis of random 
vibration using power spectrum techniques. Three integer parameters 
may be specified: 

Parameter 1: The number of discrete frequencies at which power 
spectral density values are specified. 

Parameter 2: The number of breakpoints used to define the frequency 

ranges to be used in the frequency— domain integration. 

Parameter 3: If a nonzero integer is specified, the spatial loading 

correlation function is assumed to be frequency- 

dependent, and is recomputed at each integration point 
in the frequency. Defaults to 0, i.e., S w is 

AA 

computed only once. 

Both parameters are upper bounds to the actual numbers used in the 
analysis and will be adjusted if fewer entries are specified in the 
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model data section. For further elaboration on the input required 
for pcwer spectrum analysis, please refer to the PCWER, SPECTRUM and 
BREAKPOINTS options in the model data section. See also user sub- 
routine USXX for user-defined S^. 

C.57 *PULSELQAD 

This option may be used to define conveniently a force time-history 
which may be triggered repeatedly at a number of points at specified 
times. This capability may be used to model traveling pressure wave 
phenomena. Two integer parameters are required: 

Parameter 1: The number of discrete points in time at which the 

pulse load time-history is specified. 

Parameter 2: The number of occurrences of the pulse load at 

different times and locations in the course of the 
analysis. 

Both parameters are upper bounds to the actual numbers used in the 
analysis and will be adjusted if fewer entries are specified in the 
model data section. 

C.58 *PRESSURE 

This option is used to invoke the nodal surface pressure distribu- 
tion option for meshes assembled for 2D or 3D continuum-type 
elements. These include element types 3, 7, 10, 11 and 151 through 
154. 

No parameters are needed. 

The algorithm is based on the assembly of nodal vectors correspond- 
ing to a unit pressure on all faces of the element. On a continuum 
mesh, these vectors cancel out at all interned nodes, producing non- 
zero vectors only at the external boundary nodes. The assembled 
vectors provide well-defined boundary normals, which can then be 
used for applying nodal pressures. 

C.59 *HARDENING_SD3PE 

This option sp>ecifies the number of data entries to define a 
piecewise-1 inear hardening slqpe. 

Parameter 1: The maximum number of data entries in *W3RKHARD data 

segment in the Model Data block. 

Note that *W3RKHARD and * YIELD are options available in the Model 
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Data block. *WORKHARD specifies elastic plastic constitutive be- 
havior for the entire domain. * YIELD allows users to define the 

initiation of nonlinear material response at each node. 

C.60 *LARGE_DISPLACEMENT 

This parameter card invokes the large displacement option. Note 
that the updated Lagrangian algorithm is implemented to follow the 
geometry change. 

Parameter Is Increment number at which the large displacement terms 
are turned on. 

C . 6 1 *FINITE_STRAIN 

This option invokes the finite strain terms in elastic-plastic 
computations . This parameter card cannot be specified without also 
specifying *LARGE_DISPLACEMENT card. 

Parameter 1: Increment number at which the finite strain terms are 

turned on. 

C . 62 *FOLt£WER_PORCE 

This option allows the loading condition to be updated following 
large displacement of the structure. This option cannot be used 
without also specifying the *LARGE_DISPLACEMENT card. 

Parameter 1: Increment number at which this option is turned on. 

C.63 *END 

This indicates the end of the parameter data input. 
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D. MODEL DATA 

The rrodel data follows immediately after the *END line of the parameter 
data. It consists of a number of segments. Each segment starts with a 
keyword line with at most one integer on the line. 

The number of data records following the keyword line does not have to 
be determined beforehand. The program reads data records until a new 
keyword line is encountered. 

The current version of MHOST supports the following options: 

SECTION D INDEX 


KEYWORD 

SECTION NO. 

PAGE 

*BASEXCITATIQN 

D.32 

49 

♦beamsection 

D.26 

46 

♦BODYFORCE 

D.21 

44-45 

♦BOUNDARY 

D.2 

36 

♦BREAKPOINTS 

D.39 

52-53 

♦COORDINATES 

D.8 

40 

♦DAMPING 

D.24 

45-46 

♦DASHPOTS 

D.29 

47-48 

♦DEBUG 

D. 15 

43 

♦DISTRIBUTEDLQAD 

D.19 

43-44 

♦EMATRIX 

D.34 

50 

♦DUPLICATENODE 

D.22 

45 

♦ELEMENTS 

D.l 

36 

♦END 

D. 16 

43 

♦ENDPCWER 

D.40 

53 

♦ENDPULSED3AD 

D.44 

54 

*PQLL£WER_PORCE 

D.50 

55 

♦FORCES 

D.4 

37-38 

♦HARM3NICLQAD 

D.31 

48 

♦HOLE 

D.36 

52 

♦INCLUDE 

D.45 

54 

♦INCREMENTS 

D. 10 

40-41 

♦ITERATIONS 

D.6 

38-39 

♦LAMINATES 

D.33 

49-50 

♦MASSES 

D.30 

48 

♦NODALDEFINITION 

D.43 

53-54 

♦ORIENTATION 

D.35 

50-51 

♦POWER 

D.37 

52 

♦PRESSURE 

D.49 

54-55 

♦PRINTOPTIQN 

D. 13 

41-42 

♦PROPERTIES 

D.7 

39-40 

♦PULSELOAD 

D.41 

53 

♦RESUME 

D.46 

54 
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♦SAVE D.14 43 

♦SHIFT D.27 46-47 

♦SPECTRUM D.38 52 

♦SPRINGS D.28 47 

♦STOP D.17 43 

♦STRESS D.12 41 

♦TEMPERATURES D.ll 41 

♦TIME D.23 45 

♦TIMEHISTORY D.42 53 

♦TRANSFORMATIONS D.5 38 

*T™G D.3 36-37 

♦UFXORD D.20 44 

♦WORKHARD D.18 42 
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D.l * ELEMENTS 

This data segment defines the element connectivity. When * EMBED 
option is flagged, the subelement mesh pattern can be defined by an 
additional integer entry in the data line. The number of data lines 
must be less than or preferably equal to the upper bound value given 
in the parameter data block. 

One Parameter. Containing the element type number. 

Subsequent data records each contain N+l integers, with N being the 
number of nodes per element. If the element is being further refined 
into a subelement mesh, an (N+2)nd integer entry is also required. 

Integer 1: Element number. 

Integers 2 to N+l: Node numbers of the element. 

Integer N+2 : Key for the subelement refinement. 

If the element type has more than 15 nodes continuation lines must 
be used (see B.4) . 

D.2 * BOUNDARY 

This data segment is used to impose the constraints on nodal dis- 
placements. The number of data lines must be less than or equal to 
the upper bound value given in the parameter data block. 

No Parameters. Subsequent data records each contain two (2) in- 
tegers and a real. 

Integer 1: Node number. 

Integer 2: Degree-of-freedcm. 

Real 1: Prescribed displacement value. (Default value is zero). 

If multiple specifications occur for the same degree-of-freedcm at 
the same nodal point, the last specified value is used. 

D.3 *TYXNG 

This keyword is used to set up the data segment for constraints 
between nodal displacement degrees-of-freedcm. The number of con- 
straint equations must be less than or equal to the upper bound 
value given in the parameter data block. 

Each tying constraint is specified by two data records . A typical 
form of tying constraint is: 
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D.4 


u =Sa U 
s nr in m 


where s denotes the 'slave' degree-of- freedom ('tied' degree-of- 
freedan) and m denotes the 'master' degree-of- freedom ('retained' 
degree-of freedom) . The first record specifies the nodes and de- 
grees-of-freedcm involved in the tying, the second record specifies 
the multiplication factors. 

First Record: 

Integer 1: Number of degrees -of-freedcm involved in the tying. 

Integer 2: Node number of the tied degree-of-freedan. 

Integer 3: Degree-of-freedan number of the tied degree-of- freedan. 

Node number of the first retained degree-of-freedan. 


Integer 4: 


Integer 5: Degree-of-freedan number of the first retained degree- 

of-freedan. 

Integer 6 to 2N+1: 

Node numbers /degree— of— freedan numbers of subsequent 
retained degrees-of-freedon (if N > 2). 

NOTE: Not more than 16 integers should be entered on each data line 

(see B.4) . 

Second Record: 

Real 1: Multiplication factor of the first retained degree-of- 

freedan in the tying equation. 

Real 2 to N-l: 

Multiplication factor for a subsequent retained degree- 
of-freedan in the tying equation (if N>2) . 

Real N: Prescribed value for a non-hanogeneous tying equation. 

NOTE: Not more than eight (8) reals should be entered on one data 

line (see B.4) . 

If the same degree-of-freedan appears as the tied degree-of-freedan 
in more than one tying equation, the last prescribed tying equation 
applies . Note that tying constraints take precedence over boundary 
conditions (nodal displacement constraints). 

♦FORCES 
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This keyword is used for the data segment of nodal force vectors. 
The number of data lines must be less than or equal to the upper 
bound value specified in the parameter data block. 

No Parameters. Subsequent data records each contain two (2) 
integers and one (1) real. 

Integer 1: Node number. 

Integer 2: Degree-of-freedcm. 

Real 1 : Magnitude of the applied point load. 

If multiple specifications occur for the same node and degree-of- 
freedcm the last specified magnitude is used. 

D . 5 *TRANSFORMATIONS 

This data segment specifies the coordinate transformation for the 
nodal degree-of-freedcm. Only the rotation can be specified by this 
option. 

No Parameters. Subsequent data records each contain two (2) 
integers and one (1) real 

Integer 1: Number of the node on whose coordnates the transforma- 

tion is applied. 

Integer 2: Axis number around which the coordinate system is 

rotated (if the node is tvro-dimensional, enter a zero). 

Real 1: Angle over which the coordinate system is rotated (in 

degrees) . 

If more than one coordinate transformation is applied at a node, the 
program executes all the entered transformations successively. This 
feature can be used to obtain transformations around an arbitrary 
axis . It is also possible to use tying in order to apply bou ndar y 
conditions in an arbitrary direction. 

Each subsequent transformation acts upon the last previous coor- 
dinate system defined at the node, not the original coordinate 
system used for the data input. 

Note that this option follcws the right hand coordinate convention. 
D.6 * ITERATIONS 

This data segment control the iterative solution process. When the 
BFGS option is used, the maximum number of iteration must be less 
than or equal to the number of BFGS vectors (default is 10) . When 
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the subelement iteration is invoked, convergence of the subelarent 
solution is tested with respect to the normalized relative displace- 


D.7 


(4) P reals terS ' 006 llne/ containin 9 one (1) integer and four 
Integer 1: Maximum number of iterations allowed. 


Real 1: 
Real 2: 
Real 3: 

Real 4: 


Maximum allowable relative error in the residuals. 

Maximum allowable absolute error in the residuals. 

Maximum allowable relative error in the root mean square 
of displacement (L matrix norm of displacement vector) . 

Maximum allowable relative error in the root mean square 
of strain energy associated with residual (energy semi- 
norm of the solution) . 

The default for the maximum number of iterations is 5. The default 
relative allowable error is 5.0. All other convergence 
criteria are inactive unless the tolerance is specified by the user 
with a non- zero positive number. 

co " ti J lues a^tfl all active convergence criteria are 
atisfied and the next incremental solution is entered unless the 
last increment is reached. 

When the solution fails to converge after the maximum number of 
iterations is reached, the code prints out the solution for the last 

iteration and then terminates the execution with ERROR logo on the 
output files . 

♦PROPERTIES 

This data segment is used to specify the linear elastic material 
property data at nodes. 

One (1) parameter, containing the element type. 

Subsequent data records contain two (2) integers and N reals, where 
Nis the number of properties specified for this element type in 
bection G • 

Integer 1: First node in a consecutive series of nodes (dunrny) . 

Integer 2: Last node in this consecutive series ( dunrny) . 

Reals 1 to N: 

Properties of the series of nodes. 
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NOTE: If N > 7, continuation lines must be used. (See B.4). 


The element property data must be read after the element type is 
defined in the same input data deck. 

D.8 *COORDINATES 

This data segment specifies the coordinates of nodal points in the 
mesh. The total number of entries must be less than or, preferably 
equal to, the maximum number specified in the parameter data block. 

The coordinate system in which the location of nodal points is 
defined is a global one not subjected to any transformation. 

One (1) parameter, containing the number of coordinates to be read 
in for this node. If not specified, the maximum number of coor- 
dinates for all element types as specified under C.l is used. 

Subsequent data records contain one (1) integer and N reals, where N 
is the number of coordinate directions specified or the default. 

Integer 1: Node number. 

Reals 1 to N: 

Coordinates of the node. 

NOTE: if N > 7, continuation lines must be used (see B.4) . 

For the shell element (Element Type 75), four real entries are 
required in each data line. The first three real entries are the 
global coordinates of the current nodal point with the fourth entry 
being the thickness of the shell. If the thickness is left un- 
specified, the value in the previous data line is automatically 
substituted. At least the first data line must contain the shell 
thickness entry. 

For the two-node beam element (Element Type 98), six real data 
entries are required in each data line. The first three real data 
entries are the three global node coordinates. The second three 
real data entries are the three global components of the unit vector 
defining the direction of the beam cross-section y-axis. 

D.10 * INCREMENTS 

This data segment is used to specify the naximum number of incre- 
ments. The actual number of increments are usually less than the 
value specified here. The incremental process continues until the 
maximum number is met, *STOP is encountered or no other increment 
control parameter is specified. Usually, it is not necessary to 
include this data segment in an input data file. 
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^parameters. Followed by one data line containing one (1) in- 

Integer 1: Maximum number of increments . 

The default is 30. 

D. 11 *TEMPERATURES 


The nodal point temperatures for the initial condition (zeroth 

DeratHre i^ re spec . 1 . f x ied ^ this ^ta segment. The reference tero- 
perature is zero (0) unless otherwise specified by a user 

fP°L^ SSage is issued when this model data option 
is input without the TEMPERATURE parameter card. P 

^ n th ® ® he11 elen »ent (Element Type 75), the temperature has to be 
g ven at every layer for pre-integration through the thickness. 

!^l P c ar ^ ter M .Subsequent data records contain two integers and N 
reals, where N is the number of temperatures to be specified at 
nodes. N is one for all the elements except for the shell eloient 

^To? Se^f^ ° f ^ P°**s through the t5 


D.12 


Integer Is 
Integer 2: 
Real 1 to N: 
♦STRESS 


The first node number 
The last node number 
Temperature at the nodes 


D.13 


The MHDST program allows specification of stress values at nodes. 

herv' <a 0p ^f a ^ 10n J 8 :Lnconsistent with the variational formulation and’ 
hence, it is not recommended under normal situations. The imposition 

the il£ci t pSes t ?. UndarY COndition *** lead unstable iteStion in 

^f SS ^Ponents can be prescribed at nodal 
S^enl iteration only) . The computed nodal values 

dfS? by K the prescribed values . No parameters . Subsequent 

data records each contain two (2) integers and a real. 

Integer 1: Node numbers 

Integer 2 : Stress component nu mber 

Real 1: Prescribed value 

♦PRINIOPTION 
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This option specifies the generation of lineprinter output. The 
contents of the lineprinter output is determined by subsequent data 
records. Once the print output is defined, the same output can be 
obtained in subsequent increments by entering the keyword line only. 

Upon reading of the *PRINTOPTION data, the state currently existing 
is printed at the end of the increment. A new specification can be 
entered at any stage: the old print specification is then com- 

pletely overwritten. 

The print specification consists of keywords and sequences. Ttoo 
keyword groups can be distinguished; keywords describing the quan- 
tities to be printed: 

STRESS 

STRAIN 

PLASTIC 

THERMAL 

CREEP 

INCREMENIALCISPLACEMENT 

TCTALDISPLACEMENT 

FORCE 

REACTION/RESIDUAL 

EQUIVALENT 

TEMPERATURE 

and keywords describing the finite element model 
ELEMENT 

INmaiATIONPOINT 

NODE 

LAYER 

Only the first four characters of a keyword need to be specified. A 
keyword of the first group can be followed by one to three keywords 
of the second group. A particular keyword of the second group may 
occur only once, however. The combinations must be appropriate. For 
instance, an element range should not be specified for the printing 
of displacements. 

Each keyword may be followed by 1 or 2 numbers, which specify a 
sequence. For keywords in the first group, this sequence defines 
the components to be printed. If no sequence i specified, all 
components are printed. If a keyword of the first group is not 
followed by any keywords of the second group the default printout of 
nodal values is generated. 

The number of keywords in the first group (printsets) may not be 
larger than specified in the *PRINTSET parameter segment, or, if 
absent, larger than the default of 10. 
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D.14 


D.15 


D.16 


D.17 


D.18 


D. 19 


♦SAVE 


This option controls the generation of the restart file at the end 
of increment zero. 

No parameters are specified. 

i®, a restart the new save file will overwrite the 

oia restart file. 

♦DEBUG 


This option controls the debug print. When this keyword line is 

present, all element stiffness matrices are printed in the line 
printer output file. 


The amount of information produced by this option could be 
large. Recant ended to use only for snail meshes. 


extremely 


♦END 


This indicates the end of the model data input. 
♦STOP 


With this option, program execution is terminated immediately, if 
i M the 1 u f e 5 ® intention only to check the input data, the option 
can be included before any chosen *END line. It can also be used 
for termination of plot or print specification. In that case, the 
^-is terminated immediately after generation of a plot or 


♦WORKHARD 


This segment must be entered when an elastic-plastic material is 
specified by CONSTITUTIVE in the parameter data block. 


One parameter , which specifies the number of data points used to 
specify the stress-strain law. The stress-strain relation is given 
as equivalent stress versus equivalent plastic strain, one pair of 
data points per card. Note that the first pair will specifythe 
yield stress , at zero plastic strain. A maximum of 50 data lines 
can be read in by this option. Data entered in triples as: 


REAL 1: 
REAL 2: 
REAL 3: 


equivalent 

equivalent 

equivalent 


stress 

plastic strain 

shift for kinematic and combined hardening 


♦DISTRIBUTEDLQAD 


- 43 - 



MHOST USERS' MANUAL Version 4.2 


This option defines gravity, centrifugal and traction-type dis- 
tributed loadings. The keyword has one parameter input, the element 
type. This is follcwsd by at least one data line to specify the 
loading condition. The fonret is: 

Integer 1: the first element of a series of elements subjected 

to traction loading. 

Integer 2: the last element of the series. 

Integer 3: the index for the loading type. 

REAL 1: intensity of loading type 1 for the current element. 

REAL 2: intensity of loading type 2 for the current element. 

REAL 3: intensity of loading type 3 for the current element. 

See Section G, ELEMENT LIBRARY, for details of the distributed load 
input for a given element type. 

D.20 *UFXORD 

Two integer parameters, which specify the range of node numbers for 
which user subroutine UFXORD is used to generate coordinates. 

The user subroutine must be in the load module when MHOST is ex- 
ecuted with this option. 

D.21 *BQDYFORCE 

To invoke body force loading, including both gravity and centrifugal 
loadings. One integer parameter I BODY to select an option: 

1. the gravity force 

2. the centrifugal force 

with one or two data containing two or three real parameters. 

If the parameter IB0DY=1, one data line of real numbers to specify a 
direction vector for the gravity acceleration. The data line is: 

REAL 1: x component of gravity direction vector. 

REAL 2: y component of gravity direction vector. 

REAL 3: z component of gravity direction vector. 

If the parameter IB0DY=2, two data lines of real numbers are needed 

to specify a pair of points to define the axis of rotation. Each 
line consists of a set of the coordinate data. The data line is: 

REAL 1: x coordinate of a point on the axis. 
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D.22 


D.23 


D.24 


REAL 2: 
REAL 3s 


y coordinate of a point on the axis. 
^ coordinate of a point on the axis . 


♦DUPLICATENODE 


This parameter flags the duplicate node option input, 
data records each contain 2 integers. 

Integer 1: Slave node number 


Subsequent 


Integer 2: Master node number 


♦TIME 



Integer 1; 
Integer 2: 

REAL 1: 


Maximum number of time increments. 

Flag for adaptive time step selection (non-zero value is 
required to invoke this option) . 

Size of time increment. 


REAL 2: 
♦DAMPING 


Total period of time for the analysis. 


With this option structural damping is introduced for transient 
i^^uired 1CUlati0nS ‘ 006 line containin 9 two real parameters 

REAL Is Coefficient for the mass matrix (w , see below) 

REAL 2: Coefficient for the stiffness matrix (w , see below) 

Damping is introduced to the equation of notion 
Ma+Cv+Ku=F 


where 


C - w x M + w 2 K 

Note that K is updated whenever the material tangent is recalculated 
and, therefore, C will also be updated. 

If the mode superposition option for linear dynamics analysis has 
been flagged by specifying a number in the continuation line to the 
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*MODAL option in the parameter data section, mo dal damping is as- 
sumed. In this case, instead of the format specified above, one or 
more data lines each containing two integers and one real are ex- 
pected. Data lines will specify: 

Integer 1: First mode in a consecutive series. 

Integer 2: Last mode in the same consecutive series. (Defaults to 

the first mode specified) . 

Real 1: Damping ratio for modes in the series. The resulting 

uncoupled equations of motion will be 

2 1 r n 

+ 2b w q + w q = --- /d > Ifl 
r r r^r r^r m v r* ' * 

r 

vhere b^ is the dairping ratio for the r-th node. 

D.26 *BEAMSECTIQN 

This option is used to specify one or more sets of cross-sectional 
properties for the elastic beam elements. Note that beam cross- 
sections are defined on a nodal basis and not on an element basis. 
Subsequent data records contain two integer and six real values , 
specifying: 

Integer 1: First node in a consecutive series. 

Integer 2 : last node in the same consecutive series . 

Real 1: Shear Area in the local x-direction. 

Real 2: Shear Area in the local y-direction. 

Real 3: Cross-Sectional Area. 

Real 4: Moment of Inertia about the local x-axis. 

Real 5: tfcment of Inertia about the local y-axis. 

Recil 6: Torsional Constant. 

For additional information on the definition of section properties, 
refer to the beam element description in Section G. 

D.27 *SHIFT 

Used to control the power shifts during eigenvalue extraction for 
linear dynamic analyses. At least one data line must follow, 
specifying: 


- 46 - 



MHOST USERS ' MANUAL Version 4 . 2 


Integer 1: Maximum number of eigenvalues to be extracted within the 

given frequency range. 

Real 1: lower bound for the frequency range, in cycles/time. 

Rsnl 2: Upper bound for the frequency range, in cycles/time. 

NOTE: The sum of the maxima specified in the first integer value 

for all pcwsr shifts must be less than or equal to the number 
of modes specified in the MODAL card. 

Any eigenvalues found outside the prescribed frequency ranges will 
be discarded and the program will check whether all the eigenvalues 
within each range have been found. A warning message is printed 
whenever there is the possibility of more than the specified maximum 
number of eigenvalues being present within a given frequency range. 
If more than one frequency range is prescribed, MBOST stores the new 
eigenvalues and eigenvectors immediately after all previouslv 
searched frequency ranges. 

NOTE: Specification of two overlapping frequency ranges will result 

in all modes present within the overlap zone being extracted 
and stored twice. 

D.28 *SPRINGS 

Specifies additional stiffness values at a node that are to be added 
directly into the stiffness matrix. This option may be used, for 
example, to introduce springs to the ground, or a foundation stiff- 
ness. Subsequent data lines contain: 

Integer 1: Node number for the added stiffness. 

Integer 2: Degree-of-freedom at the node. 

Integer 3: Second degree-of-freedcm at the node, used to specify 

off-diagonal terms. Defaults to the first degree-of- 
freedan specified (uncoupled springs). 

Real 1: Magnitude of the stiffness term to be added at this degree- 

of-freedom. 

NOTE: If a second degree-of-freedom is specified, it must be 

greater than or equal to the first one. That is to say, only 
the upper-triangular off-diagonal terms can be entered. 

D.29 *DASHPOTS 


Specifies additional viscous damping terms at a node, to be added 
directly into the damping matrix. Subsequent data lines contain: 
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Integer 1: Node number for the added damping. 

Integer 2: Degree-of-freedom at the node. 

Integer 3: Second degree-of-freedon at the node, used to specify 

off-diagonal terms. Defaults to the first degree— of- 
freedcm specified ( uncoupled dampers) . 

Real Is Magnitude of the viscous damping term to be added at this 
location. 

D.30 ^MASSES 

Specifies additional mass terms at a node, to be added directly into 
the mass matrix. This option may be used to include lumped masses 
not included in the finite element discretization. Subsequent data 
lines contain: 

Integer 1: Node number for the added mass. 

Integer 2: Degree-of-freedon at the node. 

Integer 3: Second degree-of-freedon at the node, used to specify 

off-diagonal terms. Defaults to the first degree-of- 
freedon specified (no coupling of mass). 

Real 1: Magnitude of the mass term to be added at this location. 

NOTE: Same as for *SPRINGS options, above. 

D.31 *HARM34ICLQAD 

This option is used to prescribe harmonic nodal force loading in the 
form of 

f = 

where (fc) is the single excitation frequency entered in the continua- 
tion line to the HARMZWICLQAD option in the parameter data block. 
Subsequent data records contain two integers and two reals, specify- 
ing: 

Integer 1: Node number for the applied load. 

Integer 2: Degree-of -freedom at the node. 

Real 1: Magnitude f of the applied load. 

Real 2: Phase angle D of the applied load. The phase angle value 

defaults to zero, which corresponds to a cos(*)t loading. 
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D.32 


D.33 


*BASEEXCITATIQN 

This option is used to prescribe harmonic base excitation in the 
torm of 


u = ue 1 ^ +A ) 

at one or more base locations, where w is the single excitation 
equency entered in the continuation line to the *BASEEXCITATION 

op ion in the parameter data section. Subsequent data records 

contain tv*D integers and tvro reals, specifying: 

Integer 1: Node number for the imposed base excitation. 

Heal 1: Amplitude u of the excitation. 

Peal 2s Phase angle D of the excitation. The phase angle value 
defaults to zero, which corresponds to cos wt excitation. 

NOTE: A base excitation can only be applied at an unconstrained 

degree -of - freedom . If the model is to be excited simul- 
taneously at all base points, it must be modeled as a free 
oody during eigenvalue extraction. 

♦LAMINATES 




This option defines material properties necessary to analy: 

^MPn^r 1 ^ 1 ^® 5 m ? del f 1 with the “HOST shell element. The 

n2^^ a E n° P V 0n in the P arameter ^a block must be input to 
properly allocate memory. No parameter. F 


C °P S1Stency of t* 1 ® t^t data stream for this option is 
^ current version of MHOST. Inconsistent specifica- 
tion of data may result in an execution error or unreliable results. 


The material description here consists of 
formats as follows: 


four data lines with 


Line 1 Integer 1: The first nodal point identification of the 

series of nodes at which the iraterial modulus is 
given by the following three data lines. 

Integer 2 : The last node of the series . 

Line 2 Reals 1 to 8: Entries of material modulus for membrane 

terms. 

Line 3 Reals 1 to 8: Entries of material modulus for membrane- 

bending coupled terms. 


- 49 - 



MHQST USERS' MANUAL Version 4.2 


Line 4 Reals 1 to 8: Entries of material modulus for bending 

terms. 

The above numeric data set must be repeated for ail nodal points. 

The entries are: 

Real 1 - local s - e component of the material moduli. 

AA XX 

2 - local s - e component of the material moduli. 

aa yy 

3 - local s^ - e^ component of the material moduli. 

4 - local s^ - e^ component of the material moduli. 

5 - local s - e component of the material moduli. 

yy *y 

6 - local s - e ccnponent of the material moduli. 

Ay xy 

7 - local s^ - e^ ccnponent of the material moduli. 

8 - local s - e ccnponent of the material moduli. 

yz yz 

where local indicates the nodal coordinate system defined at ea ch 
individual nodal point. 

D.34 *DMATRIX 

Entries to the material matrix can be read in directly by invoking 
this option for continuum elements (element types 3, 7, 10, 11, 101, 
102 and 103). No parameter. 

The material description records consist of NELSTR + 1 data lines. 
(NELSTR = number of element stress components). The format is: 

Line 1 Two integer data: 

Integer 1: The first nodal point identification number of 

nodes at which the entries of DMATRIX are given 
by the following NELSTR data lines. 

Integer 2: The last node of the series. 

Lines 2 to NELSTR NELSTR real data on each line. 

Reals 1 to NELSTR: The entries of the current row of the DMATRIX. 

D.35 ORIENTATION 
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This option is used to define the preferred orientation of the 
material axes for anisotropic materials. Subsequent data lines 
contain two integers and three reals: 

Integer 1: First node of a series of consecutive nodes. 

Integer 2: last node of a series of consecutive nodes. 

Real 1: Rotation angle , in radians 

Real 2 : Rotation angle , in radians 

Real 3: Rotation angle , in radians 

The material axes for all nodes in the series will be oriented by a 
composition of three orthogonal transformations: (1) a rotation of 

about the Z-axis, followed by (2) a rotation of f about the new y- 
axiSf and by (3) a rotation of q about the final x-axis. These 
transformations are illustrated in Figure D.l. The righthand rule 
for rotations is used throughout. 

^ ^ X, Y, Z: Global Axes 

x, y, z: Material Axes 



Z , z 

FIGURE D.l Coordinate trans formations for anisotropic 
material axes orientation 
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D.36 *HXE 


This keyword line defines an embedded hole in a globed, element. The 
* EMBED line must be included in the parameter data section. No 
parameter. 

The following data lines consist of two integer and one reeil 
variables . 


Integer 1: 


The identification number of the element in which a hole 
is embedded. 


Integer 2: The subelement type. If ein element type other than 

the global element type is specified here, the 
additional subelement type must be included in the 
Parameter Input Block. 


Real Is The dimensionless diameter of the embedded hole. The 
value is the diameter of the hole in the isoparametric 
space associated with the global element. 


The number of subdivisions in the subelement mesh for the embedded 
hole is explicitly coded in the subroutine HOLE IN. The recommended 
subelement division of two meshes along each element edge and two 
meshes in the radial direction is provided with the installation. 
This option is limited to two-dimensional analysis only. 

D.37 *POWER 


This option is used to specify the input for frequency-domain 
analysis using the power spectrum technique. It must be followed by 
the SPECTRUM and BREAKPOINTS input options and terminated with the 
ENDPCWER option . No parameters are required. 

D.38 *SPECTRUM 

This option is used to define the power spectral density function 
(PSD) to be used in the analysis. Subsequent data records contain 
two real values: 

Reeil 1: The frequency (in radians /time ) . 

Reeil 2: The PSD value at the frequency above. 

The frequencies must be input in ascending order, since linear 
interpolation will be used to evaluate PSD values at intermediate 
points. 

D.39 *BREAKPOINTS 
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This option is used to specify the breakpoints defining the fre- 
quency ranges over which the frequency-domain integration is 
performed. Subsequent data lines contain a single real value: 

Real 1: Frequency breakpoint (in radians/tirre) . 

The spectral density function for modal response will be integrated 
over the frequency range between each pair of successive breakpoints 
using 5-point Gaussian quadrature. For best results, these break- 
points should be specified close together, Tin the neighborhood of 
each undamped natural frequency of the system. 

By specifying a pair of breakpoints at^-JandCt^+J where (0 is 

a resonant frequency and d is a small deviation, the response "peak" 
for the frequency should be integrated accurately. The frequency 

breakpoints must be specified in ascending order. 

NOTE: The spatial correlation function between the Fourier 

spectra of the randan loading at two points must be specified 
via the user subroutine USXX. 

D.40 *ENDPOWER 


D.41 


This option terminates the input data block for 
snd returns control to the model data reader. 


power spectrum data 


♦PULSELQAD 


This option is used to specify the input for the pulse load defini- 
tion. It must be follc^d by the TIMEHISTORY and NODALDEFINITION 
options and terminated with the ENDFULSELQAD option. No parameters 
are required. 

D . 42 *TIMEHI STORY 


This option is used to define load time-histories. Subsequent data 
lines contain two real values: 


Real 1: Time since the start of the pulse load 

Real 2: The value of the load at the tin*> above. 

The times must be input in ascending order since linear interpola- 
tion will be used to evaluate the value of the loading at 
intermediate times. For optimal results, the pulse load time- 
history should be defined at times which are integer multiples of 
the time-step size used for the analysis. 

D . 4 3 *NOQAIDEFINITION 
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This option is used to define the locations and times at which the 
pulse load time history is to be triggered in the course of the 
analysis . Subsequent data lines contain two integers and two reals: 

Integer 1: Node number at which the pulse load will occur. 

Integer 2: Degree-of-freedcm in which the pulse load will act. 

Real 1: The time at which the pulse load time history will be 

triggered at the node. 

Real 2: A scale factor for the time history at the node. It 

defaults to the last scale factor defined. If the first 
line has no scale factor defined, it will default to a 
value of 1.0 (no scaling) . 

D.44 *ENDPULSELQAD 

This option terminates the input data block lor pulse load data and 
returns control to the model data reader. 

D.45 * INCLUDE 

This option is available in the versions of MHOST compiled under 
FORTRAN- 7 7 and allows temporary input redirection to a second input 
file. It must be followed by a line containing the file name of the 
secondary input file. Upon reaching this option, MHOST will open 
the secondary input file and continue reading input from it until a 
* RESUME card or end-of-file is found. This option should NEVER be 
used within the secondary input file, but it may be used mere than 
once in the main input file, provided that * RESUME cards are used to 
terminate all secondary input files. 

D.46 *RESUME 

This option is available in the versions of MHOST compiled under 
FORTRAN-77 and allows redirection of input back into the main input 
file (FORTRAN unit 5). This option should only be used within 
secondary input files accessed via the * INCLUDE option. 

D.49 *PRESSURE 

On meshes assembled from continuum-type elements, this option may be 
used to specify nodally defined pressures over a range of external 
boundary nodes. This option is not available for meshes comprising 
beam or shell-type elements. 

No parameters cue needed. 

Subsequent data lines contain two (2) integers and one (1) real: 
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Integer 1: The first of a series of nodes to which nodal surface 

pressures are applied. 

Integer 2 : The last node in the serie s . 

Real 1: Intensity of pressure loading, defined as positive into 

the body. Nodal pressures are specified in units of 
force per unit area on 3D and axisymmetric problems, and 
force per unit length for plane strain and plane stress 
problems . 

D . 50 *FOLLDWER_FORCE 

This keyword is used for the data segment of nodal follower force 

vectors in conjunction with the large displacement option. 


The number of data lines must be less than or equal to the upper 
bound value specified in the parameter data block. 

No Parameters. Subsequent data records each contain two (2) in- 
tegers and one (1) real. 

Integer 1: Node number. 

Integer 2: Degree-of-freedcm. 

Real 1: Magnitude of the applied point load. 

If multiple specifications occur for the same node and degree-of- 
freedom the last specified magnitude is used. 
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E. INCREMENTAL DATA INPUT 

The general ccivirents made in Section D regarding structure of the model 
data input also apply to the INCREMENTAL DATA input. All quantities entered 
in this data block, except the temperature data, are INCREMENTAL quantities, 
which are used to update existing values. If certain incremental quantities 
in this block are not specified, the program increments the quantity with 
the same value as in the preceding increment. 

The following value specified in the options are available to define 
incremental data: 


SECTION E INDEX 


KEYWORD 

♦ACCELERATION 
*AUTO INCREMENT 
♦BODYFORCE 
♦BOUNDARY 
♦DEBUG 

♦DISPLACEMENT 

♦DISTRIBUTEDLQAD 

♦END 

♦ENDINITIALCONDITICN 

*FOLLOWER_FORCE 

♦FORCES 

* INITIALCONDITIGN 

♦ITERATIONS 

♦PERIODICLOADING 

♦PRINIOPTION 

♦PROPORTIONAL 

♦SAVE 

♦STOP 

♦STRESS 

♦TEMPERATURES 

♦TIME 

♦TYING 

♦VELOCITY 


SECTION NO. 

PAGE 

E.20 

60 

E.9 

58 

E. 16 

59 

E.l 

57 

E. 11 

58 

E.18 

59 

E.4 

57 

E. 15 

59 

E.22 

60 

E.23 

61 

E.3 

57 

E.17 

59 

E.7 

57-58 

E.21 

60 

E. 12 

58-59 

E.8 

58 

E. 10 

58 

E. 14 

59 

E.6 

57 

E.5 

57 

E.13 

59 

E.2 

57 

E.19 

59-60 
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E.l * BOUNDARY 


Input is the same as in the model data block. If the specified 
boundary condition existed in a preceding increment (s) , the boundary 
condition is updated. If no boundary condition was specified in a 
preceding increment (s) , a new boundary condition is added . Boundary 
conditions can be removed by specifying the appropriate node number 
and degree-of-freedom number, with the node number preceded by a 
minus (-) sign. The boundary condition is then removed from the 
list. Removing boundary conditions creates space for additional 
boundary conditions, up to the maximum limit defined in the 
parameter data. 

E.2 *TYING 

Input is the same as in the model data block. If the tied degree- 
of-freedcm was specified in a preceding increment (s) , the tying for 
the degree -of- freedom is replaced by the newly specified tying. If 
the degree-of-freedom was not specified in the preceding 
mcrement(s) , a new TYING equation is added. 

Existing tying can be removed by specifying a tying equation with 
the relevant degree-of-freedom included as the tied degree-of- 
freedom but, without any retained degrees -of -freedom. This can be 
achieved by setting the number of degrees -of-freedcm in the tying 
input equal to one. If a tying is removed, no data record specify- 
ing coefficients is required. 

E.3 *F 1 ORCES 


Input is the same as in the model da t a block. See additional came— 
nts under E.l on updating the data. 

E.4 *DISTRIBUIEDLQAD 


Input is the same as in the model data block. See a ddi tio nal ceme- 
nts under E.l on updating the data. 

E.5 *TEMPERATURES 

I npu t is the same as in the model data block. Previously specified 
temperatures are overwritten. 

E.6 *STRESS 


Input is the same as in the model data block. See additional comme- 
nts under E.l on updating the data. 

E.7 * ITERATIONS 
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Input is the same as in the model data block. This option makes it 
possible to modify the iteration parameters. 

E.8 * PROPORTIONAL 

No parameters. This option allcws multiplication of all correctly 
applied load increments ( pressures , forces, prescribed displacement 
and non-hanogenous tying equations) with a given factor. Thermal 
loads and body forces are not affected ty this option. The multi- 
plication only applies to loads which 

1. vrere read in a previous incremental data block, and 

2. load increments which were read within the current 
block PRIOR to the PROPORTIONAL option. 

One data line containing one real. 

Real Is Multiplication factor. 

E.9 *AIHO INCREMENT 

This option invokes application of the automatic load increment 
procedure. Existing incremental loads are first modified according 
to the input entered within the same incremental data block. The 
adaptive load increment procedure, based on the arc length method, 
automatically corrects the load increment size at every iteration 
step. 

One continuation line containing two integers. 

Integer 1: Number of increments for which the specified load is 

applied. With adaptive load increments, the total load 
is equal to this number times the specified load. 

Integer 2: Maximum number of adaptive increments. If not 

specified, the increment will not be adaptive. 

E.10 *SAVE 

Input is the same as in the model data block. This option writes 
the results computed by MHOST at the end of the current increment 
onto a restart file. 

E. 11 *DEBUG 

Input is the same as in the model data block. With this option 

debug printout is generated. 

E.12 *PRINDOPTION 
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E.l 3 


E.14 


E.15 


E.16 


E.17 


E.18 


E. 19 


Input is the same as in the model data block. With this 
printer output of the MHOST results is generated. 


option line 


♦TIME 


Biter time increment and total time for either a creep or dynamic 
analysis. One data line containing two integers and two reals^ 

Integer 1: Maximum number of time increments. 

Integer 2: Flag for adaptive step selection. 

REAL 1: Size of time increment 

REAL 2: Total time period to be covered. 

♦STOP 


With this option the analysis is terminated irmediately. 
*END 


Thi s option indicates the end of the current incremental data set 

nro^?oi^ < i reine T lta - 1 data block onl y contains the *END option, a 
proportional increment of one is taken Wh<=>n ncin^ 

Ajno^ro^ENT option, one END card must be input for each of the 
♦BODYPORCE 

£^*1 ^ e / S ***** model data block. This option is used 
to define the body force load increment. 

* INITLAL2GNDITIGN 

^ tial values for transient dynamic 
computations followed by options E.18 to E. 21. This option mist be 
concluded with ENDINITIALOONDITION . No parameters. P 

♦DISPLACEMENT 

ie V 5 ® 1 tha initial nodal displacements for 

transient dynamic ccmputations . 

^ ecort ? s contain one integer to identify the nodal 
of-freedon rea ^ S/ Wlt * 1 N being the number of displacement degrees - 


♦VELOCITY 
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This option is used to prescribe the initial nodal velocities for 
transient dynamic canputations . 

Subsequent data records contain one integer to identify the nodal 
point and N reals, with N being the number of displacement degrees- 
of-freedcm. 

E.20 * ACCELERATION 

This option is used to prescribe the initial nodal accelerations for 
transient dynamic canputations . 

Subsequent data records contain one integer to identify the nodal 
point and N reals, with N being the number of displacement degrees - 
of-freedan. 

E . 2 1 *PERIODICmADING 

This option is used to prescribe periodic loading and displacement 
constraint in the form of 

f = f sinu)t 

for forces, and 

u = u sin(ajt 

for displacements. 

Subsequent data records contain three integers and two reals. These 
are: 


E.22 


Integer 1: 
Integer 2: 

Integer 3: 

Real 1: 
Real 2: 


Nodal point identification. 

Flag for the nodal force/displacement specification. 

1: nodal displacement. 

2: nodal force. 

Index for the constrained degree-of-freedan (the 
nodal force vector component) . 

Period of cyclic loading. 

Intensity of cyclic loading. 


*ENDINITIALOQNDITIQN 


This option terminates the initial condition input data block. It 
must be included if the INITIALCENDHTCN option is used. 
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E.23 *FCLLCWER_FORCE 

Input is the same as in the model data block. See additional 
nts under E.l on updating the data. 


canme- 
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F. USER SUBROUTINES 

User subroutines constitute one of the real strengths of the MHOST 
program. User subroutines allow the user to substitute his own subroutines 
for several existing routines in the program. A description of each of the 
available user subroutines is given in this manual. In addition, dis- 
cussions of special routines are also included. 

The IMPLICIT REAL*8 (A-H, 0-Z) line is required in a user subroutine 
when the code is executed on 32 bits /word machines (e.g., VAX, IBM and 
PRIME) . CRAY versions of MHOST code should not contain this declaration 
line. 

It is recommended that a user not read data cards in most of the user 
subroutines. The reason for this is that the user subroutines will be 
placed into the nonlinear analysis loop, where the user cannot know how many 
times per increment the routine will be called. 

The following user subroutines are operational in the Version 4.2 of the 
MHOST program: 


SUBROUTINE NAME PAHR 


F.6 

ANPLAS 

70-71 

F.7 

CRPLAW 

72-73 

F. 10 

UBOUN 

77 

F.5 

UCOEF 

69 

F.l 

UDERTV 

63 

F.2 

UFXORD 

65 

F.ll 

USXX 

78 

F.3 

UTEMP 

66 

F.8 

VSWELL 

74 

F.9 

WALCON 

75-76 

F.4 

WKSLP 

67-68 
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F.l SUBROUTINE UDERTV 


In the MHOST program additional element types my be defined by the 
user. This is performed by giving integers 1 through 12 in C.l and 
by coding subroutine UDERIV. This routine is called on an element 
level. The description of this routine is as follows: 


subroutine uderiv( beta, vjacob,coor,xrl,xirl, char, 

* lel , ic , j inc , j iter , nelcrd , nelnfr , nelnod , nelstr , nelint , nelchr , kw, 

* lerr , neltem, ldf , j law, igaus ) ' 

c 

C* ★ * * * * 

c 

c element derivatives routine 
c 

strain displacement matrix 
determinant 
nodal coordinates 

total displacements of the nodes of an element 
incremental displacements of the nodes of an element 
element properties 
element number 
element type 
increment number 
iteration number 

number of coordinate directions per node 
number of degrees of freedom per node 
number of nodes per element 

number of stress/ strain components per integration point 
number of integration points per element 
number of properties per element 
output device 
error flag 
not used 

product of number of nodes per element and number of 
degrees of freedom 
constitutive law type see cl.l 

if reduced itegration set to 1 , otherwise equal to nelint 
c 

implicit real *8 (a-h,o-z) 

dimension beta ( nelstr, idf , nelint) ,vjacob(nelint) , 

* coor( nelcrd, nelnod) ,xrl( nelnfr, nelnod) ,xirl( nelnfr, nelnod) , 

* char (nelchr) ’ 

c * * * * 
c 

c the user needs to define : 
c 

c beta (strain displacement matrix) for all integration points 


c 

beta 

c 

vjacob 

c 

coor 

c 

xrl 

c 

xirl 

c 

char 

c 

iel 

c 

ic 

c 

jinc 

c 

jiter 

c 

nelcrd 

c 

nelnfr 

c 

nelnod 

c 

nelstr 

c 

nelint 

c 

nelchr 

c 

kw 

c 

ierr 

c 

neltem 

c 

idf 

c 


c 

jlaw 

c 

igaus 

c 

c* * 

* * ★ * 
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c vjacob value of the jacobian at all integration points 
c 

c * * * * 

return 

end v - 
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F.2 SUBROUTINE UFXORD 


in toe MHOST program, the user may define a global coordinate system 
™?L li£ow e recta "9 ular Cartesian systan by means of user sub- 

b ®. used to "»odi«y or add to coordinate input through the 
COORDINATE option. UFXORD may also be vised. as an internal coordinate 
generator. The user must input the model definition card UFXORD. 
followed by a card giving the first and last numbers of a series of 

uraoRD for 1 LS?° RD h u iU 136 USed * The P ro 9 ram will then call 
UFXORD for each node in the series, so that the coordinates of that 

node may be modified or generated. The UFXORD option may be repeated 
as many times as necessary. 

Subroutine UFXORD is written with the following header cards: 


subroutine ufxord(xord,ncrd,n) 
c 

implicit real*8( a - h , o - z ) 
c 

dimension xord(ncrd) 
c 

c xord nodal coordinate to be defined or modified 
c nerd number of coordinate directions 

c n node number 

c 

return 

end 

where: 


NCRD 

XORD 

N 


is the number of coordinates per node. 

is the array of coordinates in the Nth node and 
will be passed in containing coordinates 
previously generated at the Nth node by 
COORDINATE, FXORD or UFXORD blocks. 

is the node number. 


One can use this routine to generate special 
(e.g. , cylindrical or spherical) or to convert 
dinate systems to a rectangular systan. 


coordinate systems 
from special coor- 
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F.3 SUBROUTINE UTEMP 

In the MHOST program the user can define temperature dependent 
material properties by adding the user subroutine UTEMP as follows: 


c . . . subroutine utemp . . . temperature-dependent material properties 
c 

subroutine utemp( value, if lag, node, temp) 
c 

c ********************************************************************** 
c 

c this user routine is used to specify material properties as 
c function of temperature for temperature-sens it ive materials 
c 

c input 

c value 

c iflag 

c 
c 
c 
c 

c node 

c temp 

c 

c output 

c 

c value 

c 
c 

c this routine is called by subroutines 'hooklw' and 'thrstn' 
c 

Q ************************************** ************ ************* ******* 

c 

implicit real* 8 (a-h,o-z) 
c 
c 

return 

end 

An example of this user subroutine is found in example problem 1.4 of 
Volume II. 


property value specified in the input deck 
flag to identify the specific property 
= 2 young's mudulus 
= 3 poisson's ratio 
= 4 coefficient of thernel expansion 
= 6 shear modulus for cubic crystal 
node number 
temperature 


property at the specified temperature as 
overwritten by the user 
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F.4 SUBROUTINE WKSLP 

This subroutine mates it possible for the user to program the yie ld 
stress and the corresponding workhardening slope dirsStly as a func 
tron of equivalent plastic strain and ta^atS.^SfiL? nledTS 
^ ° £ tte , slope of “= «3uii3ent stress !SSe£ 

S^s^cifiStiSTof JlT ^* ld £tress “T >» defied ^Lso. 

, specification of the latter is optional. If the value of the 
airrent yield is not given, the program will calculate it from the 

yi6ld V3lUe and thS WDrkhardenin g slopes defined in this 
Subroutine WKSLP is written with the following header cards: 


c 

c 

c 

c 

c 

c 

c 


subroutine wkslp(ep,sy,h,shift,tarp,intnod) 
implicit real *8 (a-h,o-z) 


ep 

sy 

h 


temp 


equivalent plastic strain 
yield stress 

slope in terms of equivalent stress / 

equivalent plastic strain 

temperature 


(user coding to be included here.) 

return 

end 


where: 


EBARPL 


HSLOPE 


SHIFT 
S YIELD 


is the current total equivalent plastic strain, 

£ P =rdE? 



is the work-hardening slope, defined as: 

d6_ 

de? 

is the amount of origin shift in the kinematic or 
combined hardening models. 

is the current yield str e ss 
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TEMP is the current temperature. 

ININOD is the current node number. 


NCTE: 


This is not 


d?T 

at 


with: 


the slope of the tensile stress-strain curve, which is 


d£ = d£? + dgP 

The user must take care to provide rate of change of stress with 
respect to plastic strain, NOT total strain. 

The user must define HSLOPE and S YIELD in this routine: EBARPL, 

TEMP and ININOD should not be changed. The call to t his subroutine 
is invoked by not specifying hardening slope at the *WORKHARD data 
section. Note that the first yield stress must be specified in the 
input data deck. The routine is called as required by the program 
during elastic -plastic calculations. .The number of times it is 
called per increment depends on the number of points going plastic, 
the nonlinearity of the work-hardening curve, and on tenperature 
dependence. 
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F.5 SUBROUTINE UCOEF 

The default for thermal expansion effects is isotropic, 
anisotropic coefficients are d e sired, the user subroutine UOOEF 
be used. UCOEF is coded and added to the system as follows: 


c-subroutine=ucoef called by 'thrstn' 

subroutine ucoef(coef,nelstr,ndi, ns hear, node, temp) 

implicit real*8( a- h , o - z ) 
c 

dimension coef(l) 
c 

c coef anisotropic coefficient of thermal expansion 
c nelstr number of strain components 

c ndi number of direct strain components 

c ns hear number of shear strain components 

c node node number 

c temp current temperature 

c 

(user coding to be included here) 

return 

end 


If 

can 
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F.6 SUBROUTINE ANPLAS 


The default in the MHOST code for the initial yield surface is 
isotropic. An anisotropic yield surface option can be obtained 
through the ANPLAS user subroutine. 


The anisotropic yield function and stress potential are assured as: 

W 6 s )2 * WlA a 3«^c-V 2+ 3a 4 T ^ * + *>6^ * & 

(R. Hill - Mathematical Theory of Plasticity, Oxford, 1950) 
where 


is the equivalent tensile yield stress for isotropic behavior. 

The user defines ratios of actual to isotropic yield (in the 
preferred orientation) in the array YFDIR for direct tension yield- 
ing, and YRSHR for yield in shear (ratio of act ual shear yield to 
/ 3, where s is the isotropic yield stress in sure shear. Then the 
a^'s above are derived as (Hill): 


a l 

1 i 1 

yrdir(2J**2 yrdir(T)**2 

1 

yrd!r(I)**5 

a 2 

1 i 1 

YRDfR(3y**5 YRDIRjiy**5 

i 

YRDIR(3y**5 

a 3 

1 + 1 
YRDIR(IJ**5 YRDiR(2j**5 

1 

YRDIR(3y**3 

a 4 

2 

YRSHR{3y**5 


a 5 

2 

YRSHR(T)**2 


a 6 

2 

YRSHR(iy**5 



subroutine anplas ( intnod, jlaw,ndi,nshear,yrdir,yrshr) 
c 

c ********************************************************************** 
c 

implicit real*8 ( a-h , o-z ) 
real *4 rwork 
c 

C ********************************************************************** 
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dimension yrdir(3) ,yrshr(3) 
c 

c n node number 

c ndi number of direct stress ccnponents 

c ns hear number of shear stress components 

c yrdir direct ratios 

c yrshr shear ratios 


(user coding to be included here) 

return 

end 
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F.7 SUBROUTINE CRPLAW 

Two types of conventional creep strain laws may be provided by the 
user. The first through user subroutine CRPLAW results in 
deviatoric creep strain; the second through' user subroutine VSWELL 
results in dilatation (volumetric) creep strain. These routines 
will be Ceil led for each nodal point. 

The program allows the user to input a general creep law through 
subroutine CRPLAW. 

The assumed form of the law is: 

e = t($, T, t, £; c , p, o/ 1 ,o^ 2 ' etc -) 
where: 

L is the equivalent creep strain rate, in uniaxial 

tension. 

6 is the current equivalent (J~) stress, normalized for 

uniaxial tension. 

T is the current total temperature, 

t is the current total time. 

c 

£, is the current total equivalent creep strain, normalized 

for uniaxial tension. 

p is the hydrostatic stress. 

are additional parameters to be employed at the user 
discretion. 

The program requires the user to program a creep law so that an 
equivalent creep strain increment is defined. An example is given 
below. 

Subroutine CRPLAW is written with the following header records: 


subroutine crplaw (eqcpi,eqprev,dum, temp, timiro, tot inc,intnod) 
c 

c ********************************************************************** 
c 

implicit real*8 ( a-h , o-z ) 
real*4 rwork 


c 
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t********* 


dimension dum(3) 


(user coding to be included here) 

return 

end 

The simplest way to define a creep strain increment fran a given 
rate law e - f(6, etc) is to multiply by t so that ^£ c =^t 
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F.8 SUBROUTINE VSWELL 

The user subroutine VSWELL allows the user to include pure swelling 
(dilatational) creep in MHOST. 

Subroutine VSWELL is written with the following header records: 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


subroutine vswell(swell,dss,terp,timinc,totinc f intnod) 
implicit real*8 (a-h,o-z) 


dimension dss( 3 ) 


swell 

dss(l) 

dss(2) 

dss(3) 

temp 

timinc 

totinc 

intnod 


is the incremental swelling strain 

is the equivalent stress 

is the dilatational stress 

is the volumentric creep strain 

is the total temperature 

is the size of the current time increment 

is the total time at the beginning of the increment 

is the node number 


(user coding to be included here) 

return 

end 
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F.9 SUBROUTINE WALCON 

The temperature dependent material constants for the Walker creep 
plasticity model are calculated in the user subroutine WALCON. An 

example used for the calculation of demonstration problen A 4 is 
given as follows: 


c 

c 

c 

c 


subroutine walcon(indep,ithem,tenpn, 

1 ee,anu,akl,ak2,anin,am,anl,an2,an3,an4,an5,an6,an7, 

2 anegO ,akind,siinf ,dnldt,dn2dt,dandt, 

3 an,alam,amu,cl,c2,c3,c4,c5) 

implicit real*8 (a-h,o-z) 

this subroutine is called by hypela to calculate all of the 
temperature dependent material constants 

except for a few changes of names, this routine is as delivered 
by b. cassenti u.t. research 

dimension tabt(6) ,eet(6) ,anut(6) ,aklt(6) ,anint<6) ,amt(6) ,anlt(6) 

d^i“^t(6r 3t(6> ' m4t(6) '^ t<6) '“ 6t(6) ' an7t(6>,ak2t(6) 

dimension akindt(6) ,siinft(6) 

data tabt/800 .dO, 1000 .d0, 1200 .d0, 1400 .d0, 1600. dO, 1800 .d0/ 

data eet/26 . d6 , 24 . d6 , 24 . d6 , 22 . 6d6 , 18 . 6d6. f 13 . 2d6/ 

data anut/0 . 322d0 , 0 . 328d0 , 0 . 334d0 , 0 . 339d0,0. 345d0, 0 . 351d0/ 

J* 3 a ^^ 0 ^ 1 ; d0 ' 75631 - d0 ' 9563l - d0 / 1 10696.d0,91505.d0,59292.d0/ 

data ak2t/0.d0,0.d0,0.d0,0.d0,0.d0,0.d0/ 

data anint/.059d0, .059d0, .079d0, .1497d0, .195d0, .223d0/ 
data amt/1 . 158d0 , 1 . 158d0 , 1 . 158d0 , 1 . 158d0 , 1 . 158d0 , 1 . 158d0/ 

data anlt/0.d0,0.d0,0.d0,0.d0,0.d0,0.d0/ 

data an2t/l.d7,1.9d7,1.5d7,2.d7,5.d6,l.d6/ 

data an3t/250 . dO , 320 . dO , 781 . 2d0 , 1178 . 6d0 , 672 . 6d0, 312 . 5d0/ 

data an4t/0.d0,0.d0,0.d0,0.d0,0.d0,0.d0/ 
data an5t/0.d0,0.d0,0.d0,0.d0,0.d0,0.d0/ 
data an6t/0.d0 / 0.d0,0.d0,0.d0,8.977d-4,2.733d-3/ 
data an7t/0.d0,0.d0,0.d0,0.d0,0.d0,0.d0/ 

data cmegOt/ 0 . dO , 0 . dO , -2000 . dO , -2000 . dO , -1434 . dO , -1200 . dO/ 

data akindt/0.d0,0.d0,0.d0,0.d0,0.d0,0.d0/ 
data Siinft/48.d3,48.d3,60.d3,l.dl0,l.dl0,l.dl0/ 

ntp=6 

ntpml=ntp-l 

tdif =tabt ( 2 ) -tabt ( 1 ) 

ll=tenpn 

12=tabt ( 1 ) -tdif 

13=tdif 


it= ( 11-12 )/13 

if ( it . It . 1 ) it=l 

if ( it . gt . ntpml ) it=ntpml 

f ac= ( tempm-tabt ( it ) ) /tdif 

ee= ( eet ( it+1 ) -eet ( it ) ) * f ac+eet ( it ) 

anu= ( anut ( it+1 ) -anut ( it ) ) * f ac+anut ( it ) 

akl= ( akl t ( it+1 ) -aklt ( it ) ) * f ac+aklt ( it ) 
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ak2= ( ak2 1 ( it+1 ) -ak2 1 ( it ) ) * f ac+ak2t ( it ) 
anin= ( anint ( it+1 ) -anint ( it ) ) *f ac+anint ( it ) 
am= ( amt ( it+1 ) -amt ( it ) ) *f ac+amt ( it ) 
anl= ( anlt ( it+1 ) -anlt ( it) ) *f ac+anlt ( it ) 
an2= ( an2 1 ( it+1 ) -an2 1 ( it ) ) * f ac+an2t ( it ) 
an3= ( an3t ( it+1 ) -an3t ( it ) ) *f ac+an3t ( it ) 
an4= ( an4t ( it+1 ) -an4t ( it) ) *f ac+an4t ( it ) 
an5= ( an5t ( it+1 ) -an5t ( it ) ) *f ac+an5t ( it ) 
an6= ( an6t ( it+1 ) -an6t ( it ) ) *f ac+an6t ( it) 
an7=(an7t ( it+1 ) -an7t ( it) ) *fac+an7t( it) 
anegO= ( anegO^ t ( it+1 ) -anegO t ( it ) ) * f ac-KanegOt ( it ) 
if (indep .eq. 0) goto 65 

akind=(akindt ( it+1 ) -akindt(it) ) *fac+akindt( it) 
siinf = ( siinf t ( it+1 ) -siinf t ( it ) ) *f ac+siinf t ( it) 

65 if (itherm .eq. 0) goto 73 

dnldt= ( anlt ( it+1 ) -anlt ( it ) ) /tdif / ( anl+1 .d-3 ) 
dn2dt= ( cin2t ( it+1 ) -an2t ( it ) ) /tdif / ( an2+l .d-3 ) 
derndt= (aneg0t( it+1 )-cmegOt( it ) ) /tdif/ (anegO+l.d-3) 

73 continue 
an=l . /anin 

edam=ee*anu/ ( ( I.d0-2.d0*anu)*(l.d0+anu) ) 

amu= ( 1 . dO-2 . dO *anu ) *alam/ ( 2 . dO *anu ) 

cl=2 . dO*arau*alam/ ( cilam+2 . d0*aniu ) 

c2=4 .dO*amu*(alam+arnu) / (alam+2 .d0*ainu) 

c3=2 .d0*amu 

c4=amu 

c5=alam 

return 

end 

Details of the Walker creep plasticity model can be found in B. N. 
Cassenti, Research and Development Program for the Development of 
Advanced Tine-Temperature Dependent Constitutes Relationships, NASA 
CR-168191 (1983), Appendix 2. 
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F.10 SUBROUTINE UBOUN 


UBOUN generates nodal displacement constraints . 
Arguments are: 


NODE (NKBC) : Array for the index of constrained nodal points 

NDOF (NKBC) : Array for the index of constrained degrees-of- 

freedcm 


DISP (NKBC) 
NKBC 


Constrained displacement (output) 
Number of displacement constraint 


RUNTIM 


Total time 


TDCNC 


Current time increment 


JINC : 

In the current 
eluded: 


Current increment number 
MHOST installation, the following example 


is 


in- 


1 


subroutine uboun(node,ndof , disp, nkbc, runtim,timinc, line) 
implicit real*8 (a-h,o-z) ,J ' 

dimension node(nkbc) ,ndof (nkbc) ,disp(nkbc> 
f ac=l . 570796327d0 
tl=runtim*f ac/15 .dO 
t2=tl+timinc*fac/15 .dO 


disp2=- . 00125*dsi_n(t2 ) 
displ=- . 00125*dsin( tl ) 
dispin=disp2 -displ 
do 1 i=l, nkbc 
nl=node(i) 
n2=ndof (i) 


if (nl.eq.2.and.n2.eq.l) disp ( i ) =dispin 

^ • e< 3* 3.and.n2.eq. 1) disp ( i ) =dispin 

continue 

return 

end 


- 77 - 



MHOST USERS' MANUAL Version 4.2 


F. 11 FUNCTION USXX 

The user-definable function USXX allows the user to specify problem- 
specific spatial correlation functions for load spectra to be used 
in frequency-dcrnain analysis of randan vibration problems . A list- 
ing of the function is provided below. 

c ... function usxx ... called from subroutine 'intsqq' 
c 

function usxx( ni, nj, coor, nnode, maxcrd, anega ) 
c 

Q *★★*★★★★★★***■*★*★★**★★★★★★*★★***★***★*★*★★****★**★★★***★★★★★*★★*★****★ 

C 

c user-defined correlation function betvreen the fourier spectra 
c of the random loading between ni and nj at frequency anega 
c 

c input 

c ni : node number i 

c nj : node number j 

c coor s node coordinates ( x, y, z... ) 

c nnode : number of nodes in the mesh 

c maxcrd : maximum number of coordinates per node 

c omega : frequency of excitation 

c 

c output 

c 

c usxx : correlation value for the load spectra 

c 

c usage 

c 

c (1) for full correlation, use function as it is, with usxx = 1 

c (2) for no correlation, remove 'c#' to get usxx = 1 if ni = nj, 

c and usxx = 0 otherwise 

c (3) for partial correlation, the appropriate correlation function 

c must be coded up by the user 

c 

Q ********************************************************************** 

c 

implicit real*8 (a-h,o-z) 
c 

dimension coor (maxcrd, nnode) 
c 

usxx = l.OdO 

c# if ( ni .ne. nj ) usxx = O.OdO 
c 

return 

end 
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G. ELEMENT LIBRARY 


cont^TnPd ^ i° n 13 a quick ^ference for the elements which are 

if current version of the MHOST program. The basic charac- 
terization of the elements is given by the parameters in the following 


TABLE G.l ELEMENT TYPE CHARACTERIZATION TABLE 


ELEMENT 


TYPE 


NELNFR NELSTR 


NELINT NELLAY NDI JLAW 

NEL£RD NELNOD NELCHR NELLY 


NET CMP 


NSHEAR 


PLANE STRESS 

3 

PLANE STRAIN 

11 

AXISYMMETRIC 

10 

3D SOLID 

7 

4 -NODE SHELL 

75 

2 -NODE BEAM 

98 

PLANE STRESS* 

101 

PLANE STRAIN* 

m 

AXISYMMETRIC* 


3D SOLID* 


PLANE STRESS 

151 

PLANE STRAIN 

152 

AXISYMMETRIC 

153 

3D SOLID 

154 


2 2 4 3 

2 2 4 4 

2 2 4 4 

3 3 8 6 

7 6 4 9 

6 6 2 6 

2 2 9 3 

2 2 9 4 

2 2 9 4 

3 3 27 6 

2 2 4 3 

2 2 4 4 

2 2 4 4 

3 3 8 6 


6 4 8 1 

6 4 8 1 

6 4 8 1 

6 8 24 1 

6 4 24 5 

6 1 12 1 

6 9 18 1 

6 9 18 1 

6 9 18 1 

6 27 B1 1 

6 4 B 1 

6 4 8 1 

6 4 8 1 

6 8 24 1 


3 2 12 

4 3 13 

4 3 14 

6 3 3 5 

5 2 3 6 

6 3 3 7 

3 2 12 

4 3 13 

4 3 14 

6 3 3 5 

3 2 12 

4 3 13 

4 3 14 

6 3 3 5 


DEFINITION OF ELEMENT PARAMETERS: 


NET CRD 

NELNFR 

NELNOD 

NELSTR 

NELCHR 

NELINT 

NELLV 

NELLAY 

NELCMP 

NDI 

NSHEAR 

JLAW 


Number of coordinates per node 
Number of degrees -of-freedcm per node 
Number of nodes per element 

Number of generalized stress and/or strain components 
Number of element material properties 
Number of 'full' integration points per elanent 
Total number of degrees-of-freedcro per elanent 
Number of integration layers through the thickness 
Number of stress and/or strain components at a point 
Number of direct stress components at a point 
Number of shear stress components at a point 
Type of constitutive model for D-matrix formulation 


Elanent types 101 through 104 are available only within the 
context of the subelement mesh refinement option. No details 
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will be given in this section. Not recctrmended for the use in 
the conventional analysis. 
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POUR-NODE PLANE STRESS ELEMENT - MHOST ELEMENT TYPE 3 


SUMMARY 

ffSUSTaf-S 

s * t V u y J T 

^ r r S 1 r a S1 eSian “»*» identified by the 

SSS^pSSSTS^T ^ at “* «- — t 

- £= ' V S' Sy J T 

and the stress vector is defined as: 

^ t0 x'6 y ' T *y1 T 

Tte^finition of stresses and strains follows standard engineering 
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II. DATA PREPARATION 

The particular data lines related to the Element Type 3 are described 

be lew. 

1. PARAMETER DATA 

♦ELEMENT Followed by an integer parameter giving the upper bound 
to the number of elements. After the keyvrord line, the 
number 3 is specified as an integer parameter. 

♦DISTRIBUTEDLQAD 

Used, without a parameter, when the distributed loading 
option is invoked. 

2. MODEL DATA 

♦DISTRIBUTEDLQAD One parameter specifying element type ( 3 ) . 

Subsequent data lines contain three integers and 
two reals. 

Integer 1: The first element of the series subjected to the 

distributed loading. 

Integer 2s The last element of the series. 

Integer 3: Pressure loading type. Value one (1) for traction 

along the element edge 1 and 2. 

REAL 1: Intensity of normal surface traction force (positive 

inward) as shown in Figure G.l. 

REAL 2: Intensity of body force loading. The total body 

force is this value times the absolute value of the 
vector defining the direction of the body force. In 
the case of centrifuged loading, the angular 
velocity (radians/time) should be input. 

♦PROPERTIES One parameter specifying element type (3). 

Subsequent data lines contain two integers and five 
reals. 

Integer 1: The first node of the series to which the current 

ma trial data applies. 

Integer 2: The last node of the series. 

REAL 1: The thickness of the material 

REAL 2: Young's modulus 


i 
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REAL 3: 

Poisson's ratio 

REAL 4: 

Thermal expansion coefficient 

REAL 5: 

Density 
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Eleoent Edge 
Traction Loading 


u 

y 


FIGURE G.l Plane Stress Elements Traction Loading Definition 
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EIGHT-NODE THREE-DIMENSIONAL SOLID ELEMENT - MHOST ELEMENT TYPE 7 
I. SUMMARY 

This element is a eight-noded isoparametric solid element with three 
eed0m Per n0dS aS shown in Fi( ^=e G.2. The degrees-of- 


U = [ U , U ,U | 
x' y z J 

in the rectangular Cartesian coordinate system and are identified bv the 
integer numbers 1 through 3 respectively. 

at nodes ^ ^ 

and the stress vector is defi ned as 
- " [ ^ x' ‘V^z' L xy '^yz'^zx ]T 

SvStiSi! i0n ° f Stresses 91x1 strains Allows standard engineering 
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II. DATA. PREPARATION 

The particular data lines for Element Type 7 are described below. 

1. PARAMETER DATA 

★ELEMENT Followed by an integer parameter giving the upper bound 
to the number of elements. After the keyword line, the 
number 7 is specified as an integer parameter. 

♦DISTRIBUTEDDOAD 

Used, without a parameter, when the distributed loading 
option is invoked. 


2. MDDEL DATA 


♦DISTRIBUTEDIOAD One parameter specifying element type ( 7 ) . 

Subsequent data lines contain three integers and 
tvro reals. 


Integer 1: The first element of the series subjected to the 

distributed loading. 


Integer 2: The last element of the series. 

Integer 3s Pressure loading type, value one (1) for traction 
along the element surface 1,2,3 and 4. 


REAL 1: Intensity of normal surface traction force (positive 

inward) as shown in Figure G.2. 

REAL 2: Intensity of body force loading. The total body 

force is this value times the absolute value of the 
vector defining the direction of the body force. In 
the case of centrifugal loading, the angular 
velocity (radians/time) should be input. 


♦PROPERTIES One parameter specifying element type (7). 

Subsequent data lines contain two integers and five 
reals. 

Integer 1: The first node of the series to which the current 

ma trial data applies. 

Integer 2: The last node of the series. 


REAL 1: Dummy entry 

REAL 2: Young's modulus 
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REAL 3: 

Poisson's ratio 

REAL 4: 

Thermal expansion coefficient 

REAL 5: 

Density 
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FOUR-NODE AXISYMMETRIC SOLID ELEMENT - MHOST ELEMENT TYPE 10 
I. SUMMARY 


This element is a four-noded isoparametric axisymmetric solid elpment- 


! = 'V u r ] T 

“* 316 identi£led b Y integer 
S ^“KTS." at «» n^es «d tte e lCTB nt integra- 


e= ie z .t t .t e ,-K a f 

and the stress vector is defined as 
- ~ *• ^z' fT zr ] T 

Snvition £ti ° n ° f Stresses and stralns standard engineering 
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II. DATA PREPARATION 

The particular data lines for Element Type 10 are described below. 

1. PARAMETER DATA 

* ELEMENT Followed by an integer parameter giving the upper bound 
to the number of elements. After the keyword line, the 
number 10 is specified as an integer parameter. 

*DISTTRIBOTEDLQAD 

Used without a parameter, when the distributed loading 
option is invoked. 

2. MDDEL DATA 

*DISTRIBUTEDLQAD One parameter specifying element type ( 10 ) . 

Subsequent data lines contain three integers and 
two reals. 

Integer 1: The first element' of the series subjected to the 

distributed loading. 

Integer 2: The last element of the series. 

Integer 3: Pressure loading type. Value one (1) for traction 

along the element edge 1 and 2. 

REAL 1: Intensity of normal surface traction force (positive 

inward) as shewn in Figure G.3. 

REAL 2: Intensity of body force loading. The total booty 

force is this value times the absolute value of the 
vector defining the direction of the booty force. In 
the case of centrifugal loading, the angular 
velocity (radians /time) should be input. 


The gravity force can only be applied in the axial (z) direction, and 
the centrifuged force can only be imposed around the z axis. 


♦PROPERTIES One parameter specifying element type (10). 

Subsequent data lines contain two integers and five 
reals. 

Integer 1: The first node of the series to which the current 

ma trial data applies. 

Integer 2: The last node of the series. 



i 
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REAL 1: 

Dunrny entry 

REAL 2: 

Young's modulus 

REAL 3: 

Poisson's ratio 

REAL 4: 

Thermal expansion coefficient 

REAL 5: 

Density 
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Body force loading possible only 
along z-axls 


FIGURE G.3 Axi symmetric Solid Element Traction Loading Definition 
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FOUR-NODE PLANE STRAIN ELEMENT - MHOST ELEMENT TYPE 11 
I. SUMMARY 

This element is a four-noded isoparametric plane strain element involv- 
ing two degrees-of-freedom per node as shown in Figure G.4. The 
degrees-of-freedan are 

T 

u = [ u , u 1 
1 x' y 1 

in the rectangular Cartesian coordinate system and are identified by the 
integer numbers 1 and 2, respectively. 

Four strain components are defined at the nodes and the element integra- 
tion points. These are 


c ) T 

and the stress vector is defined as 


6= l6 K ,6 y ,6 z .T xy ] T 

The definition of stresses and strains follows standard engineering 
convention. 
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II. DATA. PREPARATION 

The particular data lines related to the Element Type 11 are described 

be lew. 

1. PARAMETER DATA 

♦ELEMENT Followed by an integer parameter giving the upper bound 
to the number of elements. After the keyword line, the 
number 11 is specified as an integer parameter. 

♦DISTRIBUTEDDOAD 

Used, without a parameter, when the distributed loading 
option is invoked. 

2. MDDEL DATA 

♦DISTRIBUTEDLQAD One parameter specifying element type (11). 

Subsequent data lines contain three intege rs and 
two reals. 

Integer 1: The first element of the series subjected to the 

distributed loading. 

Integer 2: The last element of the series. 

Integer 3: Pressure loading type. Value one (1) for traction 

along the elanent edge 1 and 2. 

REAL 1: Intensity of normal surface traction force (positive 

inward) . 

REAL 2: Intensity of body force loading. The total body 

force is this value times the absolute value of the 
vector defining the direction of body force. In the 
case of centrifugal loading, the angular velocity 
(radians /time) should be input. 


♦PROPERTIES One parameter specifying element type (11). 

Subsequent data lines contain two integers and five 
reals. 

Integer Is The first node of the series to which the current 
ma trial data applies. 

Integer 2: The last node of the series. 

REAL Is Dumny entry 

REAL 2: Young's modulus 


- 94 - 



MHOST USERS' MANUAL 

Version 4.2 

REAL 3: 

Poisson's ratio 

REAL 4: 

Thermal expansion coefficient 

REAL 5: 

Density 
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POUR-NODE SHELL ELEMENT - MHQST ELEMENT TYPE 75 
I. SUMMARY 

This element is a four-node isoparametric shell element derived from the 
Reissner-Mindlin theory for plates and shells, involving six degrees-of- 
freedom per node. The degrees-of-freedom are 


3“ [ v V 

** 311 nMfcers 1 through 6 respectively. See 

Figure G. 5 for the gecmetrical definition. The characteristics of this 
element are similar to MARC Element 75. 

. s ^ rain components are defined at the nodes with 
respect to the local Cartesian coordinate system. These are 

- " [£ x'V Y xy 'V&V K x'Ky ^xyl 1 

and the generalized stress components associated with these entries are 

- = [N x' V V S xx- S vx' M x' V V T 


An orthogonal coordinate system is constructed at the element centroid 
and is used to define the eight strain components. 


- e y' ^xy'^xz'^yz' ] 


xy J 

The first three components are membrane strains, the fourth and fifth 
components are transverse shear strains, and the last three components 

COTp0nentS are used to £ °- «h^S~t 

< ^ Ls P 1 ^ ements ^ stresses are read in and printed out in terms of 

can P onents - ^ onl Y exception occurs when the element 

quantity output routines are invoked with the PRINTOPTION command. In 

this case information is given in the local element orthogonal coor- 
ornate system. 3 

The local coordinate system at a node is defined by calculating the 
el ^ n v t n 013 ^ vector for all the elements connected to the nodal point 
^ the J v J5 age taken to generate the normal-to-the-shell-surface 
• at node * Components of this vector are printed before 
reporting the results of the zeroth increment. The tangential coor- 
dinate system is calculated by taking the cross products of the nodal 
normai to global y-axis for the local x-axis, and the global x-axis for 
the local y-axis, respectively. Algebraic details can be found in 0. C. 
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Zienkiewicz, The Finite Element Method, 3rd Ed., McGraw-Hill, New York 
(1977) . 
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II. DATA PREPARATION 

The particular input data lines related to Element Type 75 are described 
here. For further detail, please refer to Sections B through E. 


1. PARAMETER DATA 

♦ELEMENT Follcved by an integer parameter specifying an upper 

bound for the number of elements. After the keyword 
line, the number 75 is specified as an integer 
parameter. 

♦DISTRIBUTEDLQAD Used, without a parameter, when the distributed 

loading is specified. 

♦LOUBIGNAC The third parameter associated with the stiffness 

matrix integration is not applicable to this 
element . 

2. MDDEL DATA 

♦COORDINATES Subsequent data records contain one integer to 

identify the nodal point, and three global 
coordinates for the particular node. In case of 
flat plates in the x-y plane, the thin coordinate 
can be emitted. The fourth recil data item is the 
shell thickness at the node. 

♦DISTRIBUTEDLQAD One parameter of element type (75). Subsequent data 

lines contain three integers and two real 
parameters . These are 

Integer Is First element in a consecutive series of elements. 

Integer 2: Last element in the same consecutive series of elements. 

Integer 3: 1 for the edge traction loading. 

10 for the lateral pressure loading. 

11 for the both. 

REAL Is Value for the uniform edge traction. 

REAL 2s Value for the lateral pressure loading. 

REAL 3s Intensity of the body force loading. The total body 

force is this value times the absolute value of the 
vector defining the direction of body force. In the 
case of centrifugal loading, the angular velocity 
(radians /time) should be input. 
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♦PROPERTIES 

Integer 1 
Integer 2 
REAL 1: 
REAL 2: 
REAL 3: 
REAL 4: 
REAL 5: 


With one parameter specifying element type (75) , 
subsequent data records contain two integer and three 
real parameters. 

First node in a consecutive series of nodes. 

Last node in the same consecutive series of nodes. 
Thickness of the shell. 

Young's modulus. 

Poisson's ratio. 

Coefficient of thermal expansion. 

Mass density. 


- 100 - 



MHOST USERS' MANUAL Version 4.2 


III. LINE PRINTER OUTPUT 


Stresses and strains are printed at each layer of integration through 
the thickness. Components printed are: 

* x 

2. <5 


3. -C 


xy 


4. X 


xz 


5. 


T 


yz 


The default number of layers is five, with the first and the last layers 
being the bottom and the top surfaces of the shell, respectively. The 
rest of the layers are uniformly distributed through the thickness. 
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Lateral Pressure Loading 
(Poaltlva vmlua opposite to 
eleaent z axis) 




FIGURE G.5 Shell Element Traction Loading Definitions 
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TWO-NODE TIMOSHENKO BEAM ELEMENT - MHOST ELEM ENT TYPE 98 
I. SUMMARY 

This element is a two-node linear isoparametric beam element derived 
from Timoshenko beam theory, involving six degrees -of-freedcm per node. 
The degrees -of-freedcm are 


- = 'V V V e x' e y' e z )T 

and are identified by integers from 1 through 6, respectively. See 
Figure G.6 for geometrical definition of this element. The characteris- 
tics of this element are similar to MARC Element 98. 

Six generalized strain components are defined at the nodes with respect 
to the local Cartesian coordinate system defined at the nodes . The six 
generalized strain components are 


t Y xz' V V £ z ' IC x' l V l< ^] T 

and the corresponding generalized stress components are 

i = tv v v v v v T 

At the element's single integration point, an orthogonal coordinate 
defined, having the z— axis aligned with the element axis and 
the y-axis in the plane defined by the element axis and the average of 
the normals at the two end nodes. The element stiffness matrix is 

derived in terms of the eight generalized strain components at this 
point, namely 



^xz' ^yz' ^"z' l Sc' ^y'^z^ 


with the first two components associated with transverse shear, the 
third with axial elongation, the fourth and fifth with transverse curva- 
ture, and the last one with axial twist. 
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II. DATA PREPARATION 

1. PARAMETER DATA 

♦ELEMENT Followed by an integer parameter with upper bound to 

the number of elements. After the keyword line, the 
number 98 must be specified as an the integer 
parameter. 

♦BEAMSECTION Followed by an integer parameter with upper bound to 

the number of different beam sections specified. 

*DISTRIBUTEDLOAD Used, without a parameter, when distributed loadings 

are specified. 

2. MDDEL DATA 

♦COORDINATES Subsequent data records contain one integer to 

identify the node number, followed by the three 
global coordinates of the node and three components 
of the normal vector (not necessarily a unit vector) 
used to define the plane of the local beam y-axis at 
the node. 

♦BEAMSECTION Subsequent data records contain two integers and six 

real parameters, specifying: 

Integer 1: First node in a consecutive series of nodes. 

Integer 2: Last node in the same consecutive series of nodes. 

REAL 1: Shear Area in the local x-direction. 

REAL 2: Shear Area in the local y-direction. 

REAL 3: Cross-Sectional Area. 

REAL 4: Moment of Inertia about the local x-axis. 

REAL 5: Moment of Inertia about the local y-axis. 

REAL 6: Torsional Constant. 

♦DISTRIBUTEDLQAD Followed by one integer parameter with element type 

(98) . Subsequent data lines contain three integers 
and three real parameters, specifying. 

Integer 1: First element in a consecutive series of elements. 

Integer 2 : Last element in the same consecutive series of elements . 
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Integer 3: 


Dummy parameter, any positive integer. 


REAL 1: 
REAL 2: 
REAL 3: 


Line load in the local x-direction, positive as shown in 
Figure G.6. 

Line load in the local y-direction, positive as shown in 
Figure G.6. 

Intensity of body force loading. The total body force 
is this value times the absolute value of the vector 
defining the direction of body force. In the case of 
centrifugal loading, the angular velocity (radians /time) 
should be input. 
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ELEMENT DISTRIBUTED 
LINE LOAD IN Y DIRECTION 
(POSITIVE AS SHOWN) 



ELEMENT DISTRIBUTED 
LINE LOAD IN X DIRECTION 
(POSITIVE AS SHOWN) 


FIGURE G.6 Beam Element Traction Loading Definitions 
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FOUR-NODE ASSUMED STRAIN PLANE STRESS ELEMENT - MHOST ELEMENT TYPE 151 
I. SUMMARY 

This element is a fully integrated four-node isoparametric qua- 
drilateral, based on an assumed strain field formulation. Hie element 
strains are computed in terms of a strain expansion of the form 

*1 + V - I-- B 2* 

B 1 + B 2* - K> V 
C 1 

where (x, y) is a local element cartesian coordinate system obtained by 
polar decomposition of the isoparametric mapping at the centroid of the 
element. This set assumed strain modes can represent all constant 
strain modes exactly, plus some pure bending solutions. The five 
generalized strain parameters {A^, A^, B 2 , C ]L > are related to the 

nodal displacement vector through a weak variational form. 

This element has two degrees-of-freedcm per node, namely 
T 

u = [u , u ] 

1 x' y J 

where (x, y) refer to the global cartesian coordinate system. Three 
strain components are defined at the nodes and the element integration 
points, namely 



all defined in terms of the global cartesian coordinate system. The 
three stress components are also defined in terms of the global system, 
and cue reported as 

6 = [ 6 ,6 ] T 

1 x' y xy J 

at each node and/or integration point. The definition of the stress 
ccmponents follows the standard engineering convention. 
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REAL 5: 

Mass density. 

REAL 6: 

Shear modulus (defaults to the isotropic value) . 

♦DISTRIBOTEDIDAD Followed by one parameter identifying the elenent 

type (151). Subsequent data lines contain three 
integers and two reals: 

Integer 1: 

The first of a series of elements which are subjected to 
distributed loading. 

Integer 2: 

The last element in the series. 

Integer 3: 

Surface traction type. Enter a value of one (1) if traction 
is specified on edge 1, 2 of the element. 

REAL 1: 

Intensity of normal surface traction force on edge 1, 2 of 
the element. Defined as positive inward (Figure G.7) and 
specified in units of force per unit length. 

REAL 2: 

Intensity of body force loading. For gravity loading, the 
total body force is computed by multiplying this value by 
the length of the vector defining the gravity direction. Ir 
the case of centrifugal loading, the angular velocity (in 
radians per unit time) should be specified here. 

♦PRESSURE 

Subsequent data lines contain two integers and one 
real: 

Integer 1: 

ftrst of a series of nodes to which nodal surface 
pressures are applied. 

Integer 2: 

The last node in the series. 

REAL 1: 

Intensity of pressure loading, defined as positive into the 
body and specified as a force per unit length. 
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II. DATA PREPARATION 

Specific input options which are pertinent to this element are sum- 
marized below: 

1. PARAMETER DATA 

* ELEMENT Followed by an integer parameter specifying an upper 

bound for the number of elements. The element type 
(151) is specified as an integer parameter in the 
following data line. 

♦IDUBIGNAC The third parameter, used to specify the integration 

rule for stiffness assembly, is ignored for this 
element. The element can be collapsed to form a 
triangle, provided that full integration is used for 
strain recovery and rcw-sum lumping used for 
computing the strain projection coefficients. The 
default values for this option are strongly 
recommended for all other applications. 

♦DISTRIBUTEDIOAD Used without a parameter to flag the presence of 

element distributed loads or body forces in the 
model. 

♦PRESSURE Used without a parameter to flag the presence of 

surface pressure loads specified on a nodal basis. 

2. MODEL DATA 

♦COORDINATES Subsequent data lines contain one integer 

identifying the node number, followed by two reals 
specifying the x and y coordinates for the node in 
the global cartesian coordinate system. 

♦PROPERTIES Followed by one parameter identifying the element 

type (151). Subsequent data lines contain two 
integers and up to six reals: 

Integer 1: The first of a series of nodes having the same set of 

properties . 

Integer 2: The last node in the series. 

REAL 1: A dumny parameter. 

REAL 2: Young's modulus. 

REAL 3: Poisson's ratio. 

REAL 4: Thermal expansion coefficient. 
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FIGURE G.7 Connectivity and Traction Loading for Element Type 151 
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FOUR-NODE ASSUMED STRAIN PLANE STRAIN ELEMENT - MHOST ELEMENT TYPE 152 
I. SUMMARY 


This element is a fully integrated four— node isoparametric qua- 
drilateral, based on an assumed strain field formulation. The elerent 
strains are computed in terms of a strain expansion of the form 


£ 

» t 



C 


YY 




+ - \) B 2 * 

+ B 2 * - v* A 2 y 


vdiere (x, y) is a local element cartesian coordinate system obtained by 
polar decomposition of the isoparametric mapping at the centroid of the 
element. This set of assumed strain modes can represent all constant 
strain modes exactly, plus some pure bending solutions. The five 
generalized strain parameters {A^ A^, B 2 , C ^ are related to the 

nodal displacement vector through a weak variational form, 
element has two degrees— of— freedom per node, namely 

u * <v v 

where (x, y) refer to the global cartesian coordinate system. Four 
strain components are defined at the nodes and the element integration 
points, namely 


- ^y' &z' ' 7 xy} T 

all defined in terms of the global cartesian coordinate system. The 
four stress components are also defined in terms of the global system, 
and are reported as 

- = [^x'6y< 6 z /'C xy ] T 

at each node and/or integration point. The definition of the stress 
components follows the standard engineering convention. 
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II. DATA PREPARATION 

Specific input options which are pertinent to this element are sum- 
marized belcw: 

1. PARAMETER DATA 

★ELEMENT Followed by an integer parameter specifying an upper 

bound for the number of elements. Hie element type 
(152) is specified as an integer parameter in the 
following data line. 

♦LOUBIGNAC The third parameter, used to specify the integration 

rule for stiffness assembly, is ignored for this 
element. The element can be collapsed to form a 
triangle, provided that full integration is used for 
strain recovery and row-sum lumping used for comput- 
ing the strain projection coefficients. Hie default 
values for this option are strongly recommended for 
all other applications. 

♦DISTRIBUTEDLOAD Used without a parameter to flag the presence of 

element distributed loads or body forces in the 
model. 

♦PRESSURE Used without a parameter to flag the presence of 

surface pressure loads specified on a nodal basis. 

2. MDDEL DATA 

♦COORDINATES Subsequent data lines contain one integer identify- 

ing the node number, followed by two reals 
specifying the x and y coordinates for the node in 
the global cartesian coordinate system. 

♦PROPERTIES Followed by one parameter identifying the element 

type (152). Subsequent data lines contain two 
integers and up to six reals: 


Integer 1: 

The first of a series of nodes having the same set of 
properties . 

Integer 2: 

The last node in the series. 

REAL 1: 

A dummy parameter. 

REAL 2: 

Young's modulus. 

REAL 3: 

Poisson's ratio. 

REAL 4: 

Thermal expansion coefficient. 
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REAL 5: Mass density. 

REAL 6: Shear modulus (defaults to the isotropic value). 

♦DISTRIBUTEDLQAD Followed by one parameter identifying the element 

type ( 152 ) . Subsequent data lines contain three 
integers and two reals: 


Integer 1: 
Integer 2: 


The first of a series of elements which are subjected to 
distributed loading. 

The last element in the series. 


Integer 3: 
REAL 1: 

REAL 2: 


Surface traction type. Enter a value of one (1) if traction 
is specified on edge 1, 2 of the element. 

Intensity of normal surface traction force on edge 1, 2 of 
the element. Defined as positive inward (Figure G.8) and 
specified in units of force per unit length. 

Intensity of body force loading. For gravity loading, the 
total body force is computed by multiplying this value by 
the length of the vector defining the gravity direction. In 
the case of centrifugal loading, the angular velocity (in 
radians per unit time) should be specified here. 


♦PRESSURE 


Subsequent data lines contain two integers and one 
real: 


Integer 1 : The first of a series of nodes to which nodal surface 

pressures are applied. 


Integer 2: The last node in the series. 


REAL 1: 


Intensity of pressure loading, defined as positive into the 
body and specified as a force per unit length. 
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FOUR-NODE ASSUMED STRAIN AX I SYMMETRIC ELEMENT - MHDST ELEMENT TYPE 153 
I. SUMMARY 

This element is a fully integrated four-node isoparametric qua- 
drilateral, based on an assumed strain field formulation. The element 
strains are computed in tenns of a strain expansion of the form 

= h + V - ^ ( B 1 + c i) B 2* - ^ C 2 z 
^yy = B 1 + B 2* " ^ (Aj + C j.) “ v 1 ^Y “ v> CjZ 
£§§ c i + C 2 Z ~ ^ (^1 + B i) - v* ~ ^ B 2* 

t*y " D 1 


where (x, y) is a local element orthogonal coordinate system obtained by 
polar decomposition of the isoparametric mapping at the centroid of the 
element, and z is the global symmetry axis. This set of assumed strain 
modes can represent all constant strain modes exactly, plus some pure 
bending solutions. The seven generalized strain parameters {A^, A^, B^, 

B 2' *~i' *~2' 3116 re -*- atec ^ bhe nodal displacement vector through a 
weak variational form. 


This element has two degrees -of- freedom per node, namely 
T 

u = [u , u 1 
1 z' r J 


where (z, r) refer to the global cylindrical coordinate system. Four 
strain components are defined at the nodes and the element integration 
points, namely 

- E £ z' ^ r' ^-0'^zr^ 

all defined in terms of the global cylindrical coordinate system. The 
four stress components are also defined in terms of the global system, 
and are reported as 

6= [<V6 r ,6 e ,t zr ] T 

at each node and/or integration point. The definition of the stress 
components follows the standard engineering convention. 
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II. DATA PREPARATION 

Specific input options which are pertinent to this element are sum- 
marized be lew: 

1. PARAMETER DATA 

* ELEMENT Followed by an integer parameter specifying an upper 

bound for the number of elements. The element type 
(153) is specified as an integer parameter in the 
following data line. 

♦IXXJBIGNAC The third parameter, used to specify the integration 

rule for stiffness assembly, is ignored for this 
element. The element can be collapsed to form a 
triangle, provided that full integration is used for 
strain recovery and row-sum lumping used for 
computing the strain projection coefficients. The 
default values for this option are strongly 
recarmended for all other applications. 

♦DISTRIBUTEDIOAD Used without a parameter to flag the presence of 

element distributed loads or body forces in the 
model. 

♦PRESSURE Used without a parameter to flag the presence of 

surface pressure loads specified on a no dal basis. 

2. MDDEL DATA 

♦COORDINATES Subsequent data lines contain one integer 

identifying the node number, followed by two reals 
specifying the z and r coordinates for the node in 
the global cylindrical coordinate system. 

♦PROPERTIES Followed by one parameter identifying the element 

type (153) . Subsequent data lines contain two 
integers and up to five reals : 

Integer 1: The first of a series of nodes having the same set of 

properties . 

Integer 2: The last node in the series. 

REAL 1: A dummy parameter. 

REAL 2: Young's modulus. 

REAL 3: Poisson's ratio. 

REAL 4: Thermal expansion coefficient. 
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REAL 5: Mass density. 

♦DISTRIBUTEDLQAD Followed by one parameter identifying the element 

type (153). Subsequent data lines contain three 
integers and two reals: 

Integer 1: The first of a series of elements which cure subjected to 

distributed loading. 

Integer 2 : The last element in the series . 


Integer 3: Surface traction type. Enter a value of one (1) if traction 

is specified on edge 1, 2 of the element. 

REAL 1: Intensity of normal surface traction force on edge 1, 2 of 

the element. Defined as positive inward (Figure G.9) and 
specified in units of force per unit length. 

REAL 2: Intensity of body force loading. For gravity loading, the 

total body force is computed by multiplying this value by 
the length of the vector defining the gravity direction. In 
the case of centrifugal loading, the angular velocity (in 
radians per unit time) should be specified here. The 
gravity force can only be applied in the axial direction 
(along the global z-axis), and the centrifugal force about 
the same axis. 


♦PRESSURE 


Subsequent data lines contain two integers and one 
real: 


Integer 1: 
Integer 2: 


The first of a series of nodes to which nodal surface 
pressures are applied. 

The last node in the series. 


Intensity of pressure loading, defined as positive into the 
body and specified as a force per unit area. 
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Body force loeding possible only 
along s-axls 


FIGURE G.9 Connectivity and Loading for Element Type 153 
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EIGHT NODE ASSUMED STRAIN SOLID ELEMENT - MHOST ELEMENT TYPE 154 
I. SUMMARY 


This element is a fully integrated eight-node isoparametric solid, based 
on an assumed strain field formulation. The element strains a^e com- 
puted in terms of a strain expansion of the form 


f \ + V + *3* + V* + I’v> ( B 3 + C 2 ) * - B 4 *x - -0 C 4 xy 

^yy _ B 1 + B 2* + B 3* + B 4** + II~ ( C 3 + y - v* C 4 *y - ^ 

^ **■ “ C 1 + C 2* + C 3 y + C 4^ + 1-P (A3 + B 2 ) * A 4 y» - ^ B 4 «x 

£ *y = D 1 + D 2* S* = E 1 + E 2* - F 1 + F 2 y 


where (x, y, z) is a local element cartesian coordinate system obtained 
5_ P ° 1 , ar ^ artp T/ tl0 n of ^ isoparametric mapping at the centroid of 
Thl ® set assumed strain modes can represent all constant 
train inodes exactly, plus some pure bending solutions. The fifteen 
generalized strain parameters {A 1 , A^ . . . , f 2 > are related to the nodal 

displacement vector through a weak variational form. 


This element has three degrees-of-freedon per node, namely 

u = [u , u , u 1 T 
L x' y' z J 


vrtiere ( x , y , z) refer to the global cartesian coordinate system. Six 

strain components are defined at the nodes and the element integration 
points , namely 



£y, c 2 ^ xy .v' yz ,r zx ] T 


all defined in terms of the global 
stress components are also defined 
are reported as 


cartesian coordinate system. The six 
in terms of the global system, and 


= [<VCy6 2 ,T .zj 1 


atercli node f nd/o f; integration point. The definition of the stress 
components follows the standard engineering convention. 
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II. DATA PREPARATION 

Specific input options which are pertinent to this element are 

summarized belcw: 

1. PARAMETER DATA 

♦ELEMENT Followed by an integer parameter specifying an upper 

bound for the number of elements. The element type 
(154) is specified as an integer parameter in the 
following data line. 

♦LOUBIGNAC The third parameter, used to specify the integration 

rule for stiffness assembly, is ignored for this 
element. The element can be collapsed to form a 
wedge or tetrahedron provided that full integration 
is used for strain recovery and row-sum lumping used 
for computing the strain projection coefficients. 

The default values for this option are strongly 
recommended for all other applications. 

♦DISTRIBUTEDLQAD Used without a parameter to flag the presence of 

element distributed loads or body forces in the 
model. 

♦PRESSURE Used without a parameter to flag the presence of 

surface pressure loads specified on a nodal basis. 

2. MODEL DATA 

♦COORDINATES Subsequent data lines contain one integer 

identifying the node number, followed by three 
reals specifying the x, y and z coordinates for 
the node in the global cartesian coordinate 
system. 


♦PROPERTIES 

Follcvred by one parameter identifying the element 
type ( 154 ) . Subsequent data lines contain two 
integers and up to six reals: 

Integer 1: 

The first of a series of nodes having the same set of 
properties . 

Integer 2: 

The last node in the series. 

REAL 1: 

A dummy parameter. 

REAL 2: 

Young's modulus. 

REAL 3s 

Poisson's ratio. 
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REAL 4: Thermal expansion coefficient. 

REAL 5: Mass density. 

REAL 6: Shear modulus (defaults to the isotropic value) . 

♦DISTRIBUTEDIOAD Followed by one parameter identifying the elemant 

type ( 154 ) . Subsequent data lines contain three 
integers and two reals: 

Integer 1: The first of a series of elements which cue subjected to 

distributed loading. 

Integer 2: The last element in the series. 

Integer 3: Surface traction type. Enter a value of one (1) if traction 

is specified on the 1, 2, 3 , 4 element face. 

REAL 1: Intensity of normal surface traction force on face 1, 2, 3, 

4 of the element. Defined as positive inward (Figure G.10) 
and specified in units of force per unit area. 

REAL 2: Intensity of body force loading. For gravity loading, the 

total body force is computed by multiplying this value by 
the length of the vector defining the gravity direction. In 
the case of centrifugal loading, the angular velocity (in 
radians per unit time) should be specified here. 

♦PRESSURE Subsequent data lines contain two integers and one 

real: 

Integer 1: The first of a series of nodes to which nodal surface 

pressures are applied. 

Integer 2: The last node in the series. 

REAL 1: Intensity of pressure loading, defined as positive into the 

body and specified as a force per unit area. 
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H. LINE PRINTER OUTPUT 

The following information is produced on the line printer file. 

1. MHOST logo, date and version number. On certain systems the date 
and time when the data is processed is also reported. 

2. Echo prints of data lines in parameter data section, and the core 
storage used by the data input processors . 

3. Echo prints of data lines in model data section, unless * NOECHO 
option is turned on. 

4. Summary of the parameter data and model data interpretation. 

5. The record of execution for the zeroth increment including the size 
of system matrix, the core storage used by the finite element 
solution processor followed by the condition of stiffness matrix and 
convergence of iterative procedure. 

As a measure of convergence, the program reports: 

a) Relative errors in residual. 

b) Absolution errors in residual. 

c) Root mean square error in displacement. 

d) Root mean square error in energy (which is the element-weighted- 
averaged relative error in residual ) . 

6. When convergence or the maximum number of iterations is reached the 

latest result has been printed in this file. " 

The quantities specified by the PRINTOPTICN data in the model data 
section are printed accordingly. 

a) For a vectorial quantities, its physical components are printed 
with respect to the global coordinate system except for points 
at which a coordinate transformation has been explicitly 
defined. At the node where the user has defined a coordinate 
transformation, all vector components are printed in the user 
defined coordinate system. 

b) For tensorial quantities, the von Mises intensity (the oc- 
tahedral value or the second invariant) is printed out as the 
zeroth component. 

7. When the execution fails, the program stops after printing the 

following error logo. In most cases, an explanatory error message 
is provided by MHOST. 
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CCCCCCCCCC 

CCCCCCCCCC 

CC 

CC 

cccccccc 

cccccccc 

CC 

CC 

CCCCCCCCCC 

CCCCCCCCCC 


cccccccc 

ccccccccc 

CC CC 

CC CC 

ccccccccc 

cccccccc 

CC CC 
CC CC 

CC CC 

CC CC 


cccccccc 

ccccccccc 

CC CC 

CC CC 

ccccccccc 

cccccccc 

CC CC 
CC CC 

CC CC 

CC CC 


CCCCCCCC 
CCCCCCCCCC 
CC CC 

CC CC 

CC CC 

CC CC 

CC CC 

CC CC 

CCCCCCCCCC 

cccccccc 


cccccccc 

ccccccccc 

CC CC 

CC CC 

ccccccccc 

cccccccc 

CC CC 
CC CC 

CC CC 

CC CC 


8. If a modal analysis is performed, the following information will be 

printed for each eigenvalue/eigenvector pair ; 

a) The eigenvalue number and value. 

b) The corresponding frequency in both radians per time and cycles 
per time. 

c) The corresponding eigenvector, normalized so that the L norm 
(square root of the sum of the squares of all elements) is 1.0. 

d) The generalized mass associated with this normalization of the 
eigenvector. 

9) When the subelement refinement option has been requested, the 

following information is generated and reported: 

a) At the end of the bulk data print out section the subelement 
mesh data is reported for each global element, including the 
physical and isoparametric nodal coordinates and element 
connectivities . 

b) After the global solution is reported, the same quantities are 
ob t a in ed in the subelement region and printed for each global 
element . 
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I. DOG FILE 


The progress of execution is recorded in the log file (FORTRAN Unit No. 
!>“* Prated st the terminal when the job is executed interactively on 
PRIME (Prunes), IBM (VM/CMS), and VAX (VMS) operating systems . An 
exanple is found in Section A.l of Volume II. 


For unsuccessful executions, a brief error summary is provided in the 
log f lie . Detailed error diagnostics are only provided in the line 
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J. POST PROCESSING FILE 

When *POST is entered in the parameter data block, MHOST results at the 
end of each specified increment are written on FORTRAN Unit No. 19 in 
formtted (80 bytes /record) form. No user control of the contents of 
this file is possible. 

All the information is reported at nodal points. The header records and 
the format specification for the mesh data are compatible with the MARC 
General Purpose Finite Element Program (Version K2) formatted post tape. 
The full specifications of this file can be found in MHOST Version 4.2 
Systems' Manual. 

The incremental solution data is written as both element data and no dal 
data. The element data records consist of the stress and various strain 
components along with their invariants at nodes. This data is in 6G13.6 
format. The nodal displacement and reaction force vector components are 
written in the same format as in MARC Kl. MENIAT versions number 5.0 
and higher are equipped to read the MHOST Version 4.2 post tape. 

For details of the element data specification see MHOST Systems' Manual. 
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K. MHOST INDEX 


KEYWORD 

PAGE 

KEYWORD 

PAffi 


- PARAMETER DATA - 


C.49 

♦HARMONICLOAD 

29-30 

C.l 

♦ELEMENTS 

17 

C.50 

♦BASEEXCITATION 

30 

C.2 

♦NODES 

18 

C.51 

♦COMPOSITE 

30 

C.3 

♦BOUNDARY 

19 

C.52 

♦STIFFENING 

30 

C.4 

♦TYING 

19 

C.53 

♦CENTRIFUGAIMASS 

30-31 

C.5 

♦TRANSFORMATIONS 

19 

C.54 

♦CONJUGATEGRADIENT 

31 

C.6 

♦FORCES 

19-20 

C.55 

♦NOECHO 

31 

C.7 

♦PRINTSETS 

20 

C.56 

♦PCWERSPECTRUM 

31-32 

C.8 

♦RESTART 

20 

C.57 

♦PULSELQAD 

32 

C.9 

♦TEMPERATURE 

20 

C.58 

♦PRESSURE 

32 

C.10 

♦LOUBIGNAC 

20-21 

C.59 

♦HARDENING SLOPE 

32-33 

C.ll 

♦STRESS 

21 

C.60 

♦LARGE DISPLACEMENT 

33 

C.12 

♦DYNAMIC 

21 

C.61 

♦FINITE STRAIN 

33 

C. 13 

♦OPTIMIZE 

21-22 

C.62 

*PQLLOWER_PORCE 

33 

C. 14 

♦CREEP 

22 

C.63 

♦END 

33 

C.15 

♦ANISOTROPY 

22 




C. 16 

♦MODAL 

22 


- MODEL DATA - 


C.17 

♦BUCKLE 

22-23 

D.l 

♦ELEMENTS 

36 

C.18 

♦THERMAL 

23 

D.2 

♦BOUNDARY 

36 

C.19 

♦CONSTITUTIVE 

23 

D.3 

♦TYING 

36-37 

C.20 

♦DISTRIBUTEDLOAD 

23 

D.4 

♦FORCES 

37-38 

C.21 

♦REPORT 

23-24 

D.5 

♦TRANSFORMATIONS 

38 

C.22 

♦TANGENT 

24 

D.6 

♦ITERATIONS 

38-39 

C.23 

♦UIHERM 

24 

D.7 

♦PROPERTIES 

39-40 

C.24 

♦UFORCE 

24 

D.8 

♦COORDINATES 

40 

C.25 

♦UTEMP 

24 

D.10 

♦INCREMENTS 

40-41 

C.26 

♦UCOEF 

24 

D. 11 

♦TEMPERATURES 

41 

C.27 

♦UPRESS 

24 

D. 12 

♦STRESS 

41 

C.28 

♦UHOOK 

25 

D.13 

♦PRINIOPTICN 

41-42 

C.29 

♦UDERTV 

25 

D.14 

♦SAVE 

43 

C.30 

♦SCHEME 

25 

D.15 

♦DEBUG 

43 

C.31 

♦PERIODICLQADING 

25 

D.16 

♦END 

43 

C.32 

♦DUPLICATEMODE 

25 

D.17 

♦STOP 

43 

C.33 

♦POST 

25-26 

D.18 

♦WORKHARD 

43 

C.34 

♦UBOUN 

26 

D.19 

♦DISTRIBUTEDLOAD 

43-44 

C.35 

♦BANDMATRIX 

26 

D.20 

♦UXFORD 

44 

C.36 

♦FRCNIALSOLUTION 

26 

D.21 

♦BODYFORCE 

44-45 

C.37 

♦DEFORMATION 

26-27 

D.22 

♦DUPLICATENODE 

45 

C.38 

♦EMBED 

27 

D.23 

♦TIME 

45 

C.39 

♦WKSLP 

27 

D.24 

♦DAMPING 

45-46 

C.40 

♦BEAMSECTICN 

27 

D.26 

♦BEAMSECTICN 

46 

C.41 

♦DISPLACEMENTMEIHOD 

27-28 

D.27 

♦SHIFT 

46-47 

C.42 

♦SHIFT 

28 

D.28 

♦SPRINGS 

47 

C.43 

♦BFGS 

28 

D.29 

♦DASHPOTS 

47-48 

C.44 

♦SPRINGS 

28 

D.30 

♦MASSES 

48 

C.45 

♦DASHPOTS 

28-29 

D.31 

♦HARMONICLOAD 

48 

C.46 

♦MASSES 

29 

D.32 

♦BASEXCITATION 

49 

C.47 

♦LINESEARCH 

29 

D.33 

♦LAMINATES 

49-50 

C.48 

♦SECANINEWION 

29 

D.34 

♦DMATRIX 

50 
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KEYWORD PACT! 

D.35 *ORIENIATION 50-51 

D.36 *HXE 52 

D.37 *POWER 52 

D.38 *SPECTRUM 52 

D.39 * BREAKPOINTS 52-53 

D.40 *ENDPOWER 53 

D.41 *PULSEDOAD 53 

D.42 *TIMEHI STORY 53 

D.43 *NODALDEFINITION 53-54 

D.44 *ENDPULSELQAD 54 

D.45 * INCLUDE 54 

D.46 *RESUME 54 

D.49 *PRESSURE 54-55 

D. 50 FQLLCWER_PORCE 55 

- INCREMENTAL DATA - 

E. l *BOUNDARY 57 

E.2 *TYING 57 

E.3 *PORCES 57 

E.4 *DISTRIBOTEDLQAD 57 

E.5 ^TEMPERATURES 57 

E .6 *STRESS 57 

E.7 * ITERATIONS 57-58 

E .8 *PROPORTIQNAL 58 

E.9 * AUTO INCREMENT 58 

E.10 *SAVE 58 

E. 11 *DEBUG 58 

E.12 *PRINrOPTION 58-59 

E.13 *TIME 59 

E.14 *STOP 59 

E.15 *EM) 59 

E.16 *BODYFORCE 59 

E.17 * INITIALCQNDITION 59 

E.18 *DISPLACEMENT 59 

E. 19 * VELOCITY 59_60 

E.20 * ACCELERATION 60 

E.21 *PERIODICLQADING 60 

E.22 *ENDINITIALCCNDITICN 60 

E.23 *POLLOWER_FORCE 61 
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